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We revisit the Twin Higgs scenario as a “dark” solution to the little hierarchy problem, 
identify the structure of a minimal model and its viable parameter space, and analyze its 
collider implications. In this model, dark naturalness generally leads to Hidden Valley 
phenomenology. The twin particles, including the top partner, are all Standard-Model- 
neutral, but naturalness favors the existence of twin strong interactions - an asymptotically- 
free force that confines not far above the Standard Model QCD scale - and a Higgs portal 
interaction. We show that, taken together, these typically give rise to exotic decays of 
the Higgs to twin hadrons. Across a substantial portion of the parameter space, certain 
twin hadrons have visible and often displaced decays, providing a potentially striking LHC 
signature. We briefly discuss appropriate experimental search strategies. 
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I. INTRODUCTION 

The principle of Naturalness, the notion that the weak scale should be insensitive to 
quantum effects from physics at much higher mass scales, necessitates new TeV-scale physics 
beyond the Standard Model (SM). It has motivated a broad program of searches at the LHC, 
as well as lower-energy precision/flavor/CP experiments. The absence thus far of any signals 
in these experiments has disfavored the most popular scenarios, including supersymmetry 
(SUSY), Composite Higgs and Extra Dimensions, unless their mass scales are raised above 
natural expectations, leading to sub-percent level electroweak fine-tuning in complete mod¬ 
els. The puzzle over why Nature should be tuned at this level has been dubbed the “Little 
Hierarchy Problem,” and it has led to a major rethinking of naturalness and its implications 
for experiments. 

In a bottom-up approach to this problem, one may take a relatively agnostic view of very 
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high energy physics, and focus instead on naturalness of just the “little hierarchy”, from an 
experimental cutoff of about 5-10 TeV or so down to the weak scale. Unlike naturalness 
considerations involving extremely high scales, such as the Planck scale, which are tied 
to multi-loop virtual effects on the Higgs sector of all SM particles, the bottom-up little 
hierarchy problem is simpler, relating predominantly to one-loop effects of just the heaviest 
SM particles, i.e. those coupling most strongly to the Higgs. In Little Higgs, (Natural) 
SUSY, and extra-dimensional models of gauge-Higgs unification (including warped models 
that are dual to Higgs Compositeness via the AdS/CFT correspondence), large one-loop 
radiative corrections to the Higgs from the heaviest SM particles cancel algebraically against 
those of new symmetry “partners” of these heavy particles, thereby ensuring stability of the 
little hierarchy. 

The most significant of these corrections is associated to the top quark. Naturalness 
requires that this must be substantially canceled by a corresponding correction from the 
top’s partner(s), to which it is related by an ordinary global symmetry or supersymmetry. 
This requirement can only be fulfilled naturally if the associated partner has a mass scale 
~ 500 GeV, easily within the kinematic reach of the LHC. Such particles also have significant 
LHC production cross-sections, since a top partner generally carries the same color charge 
as the top quark. Through this logic, the search for top partners in the above incarnations, 
under a variety of assumptions about possible decay modes, has become a central pursuit 
of the LHC. 

But must the top partner be colored? The answer is obviously critical to experimental ex¬ 
ploration. Naively the answer is yes, because the algebraic cancellations depend at one-loop 
on the top Yukawa coupling to the Higgs, and this coupling is itself corrected at one higher 
loop by QCD. In a rare exception to the “one-loop rule”, such two-loop radiative corrections 
to the Higgs are still quantitatively important for the little hierarchy problem because of 
the strength of QCD. It would seem then that the top partner should also be colored so as 
to parametrically “know” about this QCD effect in symmetry-enforced cancellations. 

Yet, remarkably, there do exist solutions to the little hierarchy problem, “Twin Higgs” 
being the first and prime example, in which the top partners are uncolored . 1 Here 

1 Another known solution to the little hierarchy problem which involves uncolored top partners is “Folded 
Supersymmetry” [4]. 
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the cancellation among radiative corrections is enforced by a discrete Z 2 symmetry that 
exchanges SM particles with new “twin” states. One way to assure naturalness cancellation 
of QCD two-loop effects is to have the twin symmetry exchange SM color and its gluons 
with a distinct twin color gauge group and its twin gluons, which couple to and correct the 
twin top partner just as QCD does the top quark. We will focus on theories of this type, 
specifically ones in which all twin particles are “dark”, with 110 SM quantum numbers. 

Colorless twin tops are vastly more difficult to produce at the LHC than top partners 
of more popular theories; indeed this is true for all twin sector particles. Twin particle 
production can only proceed through a “Higgs portal”, a modest mixing between the SM 
and twin Higgs sectors that is a necessary consequence of the twin Higgs mechanism for 
addressing the little hierarchy. Not only is the production rate small, the hidden particles 
barely interact with ordinary matter, and (at least naively) one would expect they escape the 
detectors unobserved. The resulting missing energy signature with a very low cross-section 
would pose great difficulties at the LHC. 

How else can the twin sector be detected? Higgs mixing and virtual twin top loops 
can also subtly affect the SM-like Higgs, making precision tests of its properties extremely 
important. But with expected LHC precision, the visibility of the twin sector in this manner 
at the LHC is limited [5] . The Higgs of the twin sector can also potentially be produced, but 
it may be too wide to observe as a resonance. At best, it is heavy and has a low cross-section, 
and is far from being excluded or discovered. For these reasons, Twin Higgs remains a viable 
resolution of the little hierarchy problem, a well-hidden outpost of naturalness. 

In this paper, we re-examine the Twin Higgs scenario as an important and distinctive case 
study in “dark” naturalness, and we identify exciting new experimental opportunities for its 
discovery. We develop minimal Twin Higgs models addressing the little hierarchy problem, 
roughly paralleling the way in which “Natural SUSY” [MI] has emerged as a minimal 
phenomenological approach to the little hierarchy problem in the SUSY paradigm. I 11 both 
cases, minimalism can be viewed as a tentative organizing principle for doing phenomenology, 
starting with searches for the minimal natural spectrum and then “radiating outwards” to 
include searches complicated by non-minimal states. Also in parallel to SUSY, we find that 
the two-loop relevance of QCD interactions to the little hierarchy problem leads to some of 
the most promising experimental signals. 

In SUSY, the symmetry cancellation at two loops requires the presence of a gluon-partner, 
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FIG. 1: Example of a Twin Higgs collider event. The SM-like Higgs decays through a loop of 
the twin tops into a pair of twin gluons, which subsequently hadronize to produce various twin 
glueballs. While some glueballs are stable at the collider scale, Go+ decay to Standard Model 
particles is sufficiently fast to give LHC-observable effects, including possible displaced vertices. 
The hgg coupling, indicated by a black dot, is generated by small mixing of the Higgs and the twin 
Higgs. 

the gluino. With large color charge and spin, the glnino is phenomenologically striking over 
much of motivated parameter space, almost independent of its decay modes [T2] I14| . In Twin 
Higgs models, the analogous two-loop role is played by twin gluons, which can again give rise 
to striking signatures over a large part of parameter space, not because of large cross-sections 
but because they, along with any light twin matter, are confined into bound states: twin 
hadrons. Together with the Higgs portal connecting the SM and twin sectors, the presence 
of metastable hadrons sets up classic “confining Hidden Valley” phenomenology PH, 
now in a plot directly linked to naturalness. 

A prototypical new physics event is illustrated in Fig. [lj The scalar line represents the 
recently discovered 125 GeV Higgs scalar. This particle is primarily the SM Higgs with 
a small admixture of twin Higgs; it is readily produced by gluon fusion. But because of 
its twin Higgs content, it has at least one exotic decay mode into twin gluons, induced 
by twin top loops, with a branching fraction of order 0.1%. The twin gluons ultimately 
hadronize into twin glueballs, which have mass in the ~ 1 — 100 GeV range within the 
minimal model. While most twin glueballs have very long lifetimes and escape the detector 
as missing energy, the lightest 0 ++ twin glueball has the right quantum numbers to mix with 
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the SM Higgs, allowing it to decay back to the SM on detector timescales. The first excited 
0 ++ state also may have this property. This type of effect was first studied, in the context 
of a Hidden Valley/quirk model [121 T7, EE EEj by Jnknevich [23]. 

If the lightest 0 ++ glueball, which we call G' 0 +, has a high mass, then its decay is prompt, 
and its production rate in Higgs decays may be too rare for it to be observed at the LHC. 
However, these non-SM decays of the Higgs would be discoverable at e + e _ machines, pro¬ 
viding strong motivation for 1LC, FCC-ee, or CEPC. 

For lower mass (typically below ~ 40 GeV in the minimal model we describe) the decay of 
the G 0 + may be macroscopically displaced from the interaction point. Such displaced decays 
are a striking signature, spectacular enough to compensate for the relatively low production 
rate, and represent an excellent opportunity for the LHC. 

Moreover, the signal may be enhanced and/or enriched if there is a sufficiently light twin 
bottom quark. Depending on this quark’s mass, twin bottom production may lead to a larger 
overall twin hadron production rate; the resulting final states may include twin glueballs, 
twin bottomonium, or both. In some portions of parameter space, all twin bottomonium 
states are invisible. In others the lightest 0 ++ state (which we will call y , since it is analogous 
to the Xo states of SM quarkonium) can decay visibly by mixing with the Higgs just like 
the G*o+, though possibly with a small branching fraction. There is also the potential for 
displaced vertices from this state, though our calculations of lifetimes and branching fractions 
suffer from large uncertainties. 

With a branching fraction of the Higgs to twin hadrons of order 0.1% or greater, our 
minimal Twin Higgs model should motivate further experimental searches for this signal 
of hidden naturalness. A branching fraction of 0.1% represents ~ 500 events produced at 
ATLAS and at CMS during Run I, and we can expect many more at Run II. But triggering 
inefficiencies threaten this signal. Despite the low backgrounds for highly displaced vertices, 
triggering on such events can be a significant challenge US MM\. There is urgency as 
we approach Run II, since it is not trivial to design triggers that will efficiently capture 
all variants of the displaced vertex signature arising within the parameter space of the 
model. In particular, since twin hadron lifetimes are a strong function of their mass, and 
since displaced decays in different subdetectors require quite different trigger and analysis 
methods, a considerable variety of approaches will be required for efficient coverage of the 
parameter space of the model. Further complicating the matter is that precise theoretical 


7 


predictions of twin glueball and especially twin qnarkoninm production in Higgs decays is 
extremely difficult. Hadronization in the twin sector is a complex and poorly understood 
process, and considerably more investigation will be needed before predictions of glueball and 
quarkonium multiplicity and kinematic distributions could be possible. Thus this solution to 
the naturalness problem requires further work on both experimental and theoretical fronts. 

We should note that every element that goes into this story has appeared previously in 
the literature. Hidden glueballs appear in quirk models and other Hidden Valley models m 
[22]; they specifically arise in Folded Supersymmetry [3] and the Quirky Little Higgs 
, which attempt to address the hierarchy problem; the mixing of the G 0 + with the Higgs 
to generate a lifetime that is short enough to be observed at the LHC but long enough 
to often be displaced has appeared in a study [23J of a Hidden Valley model with quirks 
m- However we believe that this is the first time these elements have all been assembled 
together, giving the striking observable signal of exotic Higgs decays, possibly to long-lived 
particles, as a sign of hidden naturalness. 

Returning from the phenomenology to broader considerations, we note that the Twin 
Higgs shares with Little Higgs and Composite Higgs models the realization of the Higgs 
as a pseudo-Goldstone boson, while representing a significant break with earlier thinking 
and phenomenology in having a colorless top partner. In the far UV, it is possible that the 
Twin Higgs structure might match on to more conventional SUSY [33H35] or Composite 
dynamics SUEZ], or perhaps to something quite novel [351 SS]- We hope that studying the 
UV matching is facilitated by our bottom-up exploration of minimal Twin Higgs structure 
and phenomenology. We also hope that this work broadens our perspective on naturalness, 
motivates more careful investigations of other “darkly” natural mechanisms, such as “Folded 
Supersymmetry” [3], and perhaps inspires entirely new mechanisms. Ultimately, we hope to 
broaden the scope of experimental strategies motivated by naturalness, perhaps even leading 
to a discovery unanticipated by theory. 

The paper is organized as follows. Section |TT] reviews the basic Twin Higgs structure 


and develops the minimal model. We show in Section III that the electroweak tuning in 


the minimal Twin Higgs model is very mild. Section IV shows quantitatively that the twin 


sector likely contains twin QCD in order to maintain naturalness, and its confinement and 
hadrons are discussed. Section [V] derives some properties of the resulting twin hadrons, and 
discusses their production and decays. Some subtle points on hadron decays and hadron 




production are left for Appendices [A] and [B] Section |VT| synthesizes the earlier considerations 
and discusses LHC phenomenological implications. Possible experimental strategies for long- 
lived particle searches are briefly considered in Section VIA and in more detail in Appendix 
0 Precision Higgs measurements are considered in Section |VID| , and precision electroweak 
constraints are discussed in Appendix[DJ Our conclusions appear in Section VII In Appendix 
[E[ the phenomenological effect of gauging twin hypercharge, as a non-minimal extension of 
the model, is considered. 


II. THE MINIMAL OR “FRATERNAL” TWIN HIGGS 

In this section we construct the minimal Twin Higgs model, starting by reviewing the 
basic symmetries and Higgs structure and then justifying each addition to the twin sector 
based on the need to maintain naturalness and internal consistency. Because the minimal 
model does not duplicate all SM states in the twin sector (in contrast to the original mirror 
Twin Higgs model and its descendants im in which the twin sector and its couplings are 
an exact copy of the SM), we will refer to this construction as the “Fraternal Twin Higgs”. 
The model is summarized in subsection III Cl 


A. The Central Mechanism 

At its heart, the Twin Higgs mechanism involves realizing the SM-like Higgs as a pseudo- 
Goldstone boson of an approximate global symmetry, namely 517(4). An SU (4)-fundamental 
complex scalar "H with potential, 

V = A(|"H| 2 — / 2 /2) 2 , (1) 

acquires a vacuum expectation value (T-L) = f / \/2, breaking 517(4) —> 517(3) and giving rise 
to seven Goldstone bosons, of which one is ultimately identified as the SM-like Higgs scalar. 
The 517(4) is explicitly broken by the gauge and Yukawa couplings of the Standard Model. 
Without additional recourse, this explicit breaking would lead to quadratic sensitivity to 
higher scales at one loop. 

In the context of conventional global symmetry protection, this explicit breaking and UV 
sensitivity can be ameliorated by extending Standard Model gauge bosons and fermions into 
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representations of the full global symmetry, at the price of introducing additional partner 
states charged under the SM. In the Twin Higgs, the key insight is that SM states need not 
be extended into full representations of the global symmetry, and instead are merely related 
to partner states by a Z 2 exchange symmetry. This Z 2 is then promoted, at the level of 
quadratically divergent radiative corrections, to an accidental SU (4) symmetry. The partner 
particles are no longer related to SM states by a continuous symmetry, and so need not carry 
SM gauge quantum numbers. 

As in any global symmetry mechanism that stabilizes the weak scale, the Twin Higgs 
does not address the big hierarchy problem all the way to the Planck scale, but merely a 
little hierarchy up to a cutoff A, which is typically ~ 5 — 10 TeV. This also roughly matches 
the maximal reach of the LHC. Above A we imagine that one of the canonical solutions to 
the “big” hierarchy problem, such as supersymmetry or compositeness, kicks in to provide 
protection against yet higher scales. 

We can embed the SM Higgs into ([TJ) by decomposing TL = ( A,B ) into two doublets 
A and B. We identify A with the SM Higgs doublet of gauged electroweak SU(2) L (and 
charged under U(l)y of course), while B is a doublet of a different SU(2) symmetry, which 
we will call twin SU( 2). At this stage, the twin SU(2) group can be either global or gauged. 
The Z 2 symmetry acts by exchanging A -H- B. By far the largest source of explicit breaking 
from the Standard Model will be the top Yukawa, so to see the magic of the twin mechanism 
let us also introduce twin top multiplets: three species of the twin fermion Q a transforming 
as doublets of the twin SU (2) as well as twin right-handed tops u a , where a = 1 • • • 3. The 
Z 2 symmetry implies a twin top Yukawa coupling 

C D y t BQ a u a , (2) 


where we expect i/t ~ yt for the twin mechanism to be effective. Note that the a index 
implies at least a global SU{3) symmetry acting on fermions Q and it. 

The one-loop radiative potential for the A and B multiplets coming from loops of top 
and twin top quarks takes the form 

16n 2 V {1 - loop) = -6y?A 2 \A\ 2 - 6$A 2 \B\ 2 + 3^ 4 |A| 4 log ( a2 /% 2 |a| 2 ) + 3% 4 |H[ 4 log ( a2 /% 2 |b| 2 ) (3) 


for a uniform cutoff A. 2 Notice that if y t 


y t then the terms quadratic in the cutoff arrange 


2 


One should interpret this cutoff as merely a proxy for physical effects; in a realistic UV completion A will 
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themselves into an 577 (4) invariant oc |'H| 2 , while the radiative quartics explicitly break the 
577(4) but preserve the Z 2 . This is the magic of the twin mechanism: if the SM-like Higgs 
can be identified with a pseudo-Goldstone of the spontaneously broken 5f/(4), its mass will 
be insensitive to the SU (4)-symmetric quadratic divergences at one loop (or more) provided 
the Z 2 relates y t and y t . 

A second central issue is vacuum alignment. The vacuum expectation value (vev) of 
H, f/V2, which breaks the approximate global SU(4) symmetry down to SU( 3), can be 
decomposed, 

f = v\ + vl . (4) 

In general this vev breaks both SU (2)^x7(1 )y and twin SU (2), such that most components 
of "H are eaten if twin SU (2) is gauged. The remaining physical scalar states consist of linear 
combinations of the radial mode and a single uneaten Goldstone boson. We will identify 
the observed SM-like Higgs with the uneaten Goldstone boson of SU (4) / SU (3). The mass 
of the radial mode (corresponding to \T-L\ 2 fluctuations) is y/2Xf. (Alternately, we could 
work purely in terms of the non-linear sigma model of SU(4)/SU(3), in which case there is 
no perturbative radial mode, and A is identified with the cutoff of the non-renormalizable 
theory. We do not follow this approach here.) To obtain a realistic vacuum, we must break 
Z 2 to a small extent. In an exact Z 2 model the potential favors va — vr and the pseudo- 
Goldstone Higgs is an equal mixture of A and B , where recall only A carries SM Higgs 
quantum numbers. Since we have observed a SM-like Higgs experimentally, we will need 
perturbations stabilizing v\ <C so that the pseudo-Goldstone Higgs is primarily aligned 
with A. 

B. Minimal Particle Content 

Thus far we have seen how the essential structure of the twin mechanism operates at 
the level of the Higgs and the largest source of explicit SU (4) breaking, the top Yukawa 
coupling, and understood the new twin states thereby required. We continue this process of 
deducing the minimal ingredients for a realistic and natural Twin Higgs model. 

We begin with the top yukawa itself. We have seen schematically that the Higgs is 

be replaced by Z 2 -symmetric physical thresholds. 
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FIG. 2: Cancellation of the top divergence in the Twin Higgs model. The effective vertex in the 
second diagram arises upon integrating out the heavy radial mode. 


protected against large cutoff sensitivity from top loops provided a twin top with y t = y t . 
But how much can the top Yukawa Z 2 be relaxed while preserving the naturalness of the 
weak scale? When the coupling (J2]) is introduced, by Eq. (J3]) the physical mass of the 
pseudo-Goldstone boson Higgs gets a quadratically divergent radiative correction at the 
scale A, 

6m 1 ~ ^2 (y't ~ yt) ■ ( 5 ) 


This precisely cancels out when the Z 2 symmetry is exact. We can picture this in Fig. [2j 
where the pseudo-Goldstone Higgs acquires an effective coupling to the twin top upon inte¬ 
grating out the heavy radial Higgs mode. The cancellation in Fig. [2] is very similar to that 
in Little Higgs theories [40l 142] , 3 with the difference that the top partner is uncolored. But 
without exact Z 2 symmetry, the naturalness demand that these corrections are not much 
larger than the observed SM-likc Higgs mass-squared of (125 GeV) 2 , translates into 


fe(A)-yt(A) 
yt( A) 


<0.01, 


( 6 ) 


for A ~ 5 TeV. 

Ordered by the size of tree-level couplings to the Higgs doublet, the next ingredients to 
consider are the potential twin sector equivalents of SU(2) L x U(l)y gauge bosons. The 
contribution to the Higgs mass-squared from SU(2) L boson loops in the SM is ~ (400 GeV) 2 
for a cutoff A ~ 5 TeV — still a significant source of electroweak destabilization, although 
subdominant to top loops. This suggests gauging the twin SU (2) global symmetry acting on 
B and Q. Introducing twin weak gauge bosons with coupling g 2 translates to the quadratic 


3 See also [43] for review and references therein. 
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cutoff sensitivity in the pseudo-Goldstone Higgs, 


5ml 


9A 2 

32^ 


m a) - 52(A)), 


(7) 


in analogy with Eq. (J5]) . Demanding this not significantly exceed the observed Higgs mass- 

< 0.1. Note that with this gauging of twin SU(2 ) all Goldstone 


squared implies 


<72 ( A )~32(A) 


32(A) 


bosons of SU(4:) breaking, except for the SM-like Higgs itself, are now longitudinal weak 
bosons of the visible and twin sectors. 

In contrast, the contribution to m 2 from U( l)y loops in the SM is comparable to rrij, 
for a cutoff A ~ 5 TeV, and thus already consistent with naturalness. Thus naturalness 
does not require twin hypercharge, although it was included in the original Twin Higgs [1]. 
This is analogous to the statement that in natural supersymmetry there is no need for the 
Bino to be light; its presence in the low-energy spectrum is non-minimal from the bottom- 
up point of view. Given that our principle in this paper is to seek the most economical 
version of the twin Higgs that is consistent with the naturalness of the little hierarchy, we 
do not include twin hypercharge in the minimal twin Higgs model, assuming instead that 
it was never gauged or that it was broken at or around the scale A. However, for the sake 
of completeness, we will briefly discuss the significant phenomenological consequences of a 
light twin hypercharge boson in Appendix |E| 

Next we turn to the twin analogue of QCD. Of course the Higgs does not couple to SU(3) 
at tree level, but rather at one loop via its coupling to the top quark. This nonetheless leads 
to sizable two-loop corrections to the Higgs mass from physics around the cutoff. As we will 
discuss in detail in Section IV A the contribution to the Higgs mass-squared from two-loop 
y 2 g\ corrections in the SM is at least ~ (350 GeV ) 2 for a cutoff A ~ 5 TeV, putting QCD 
on similar footing as the weak gauge group. Gauging the twin SU(3) global symmetry with 
coupling g 3 gives quadratic cutoff sensitivity in the pseudo-Goldstone Higgs 

3</?A 2 


Srrif ^ 


47T 4 


(<?3 


gl) 


( 8 ) 


This is a key observation that will drive the phenomenology of a viable Twin Higgs model: 
naturalness and minimality favor a confining gauge symmetry in the hidden sector, “twin 
glue”. This twin glue has a coupling close to the QCD coupling — we will see how close in 
and therefore it confines at a scale A 3 which is logarithmically close to the 


Section IV A 


SM QCD confinement scale. 
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Once the twin SU( 2 ) and SU( 3) are gauged, we should include a variety of twin fermions 
in addition to the twin top quark to cancel anomalies. To render the twin SU( 3) anomaly- 
free we should include a twin RH bottom quark 6 ; symmetries then admit the hidden sector 
bottom Yukawa coupling 

C = y b BQb . (9) 

Unlike the case of the top sector, we should not necessarily demand that y b = y b since the 
bottom Yukawa has a much weaker effect on the SM Higgs mass at one loop. At this stage it 
suffices that y b <C yt in order to avoid creating a hierarchy problem from the bottom sector. 

Similarly, canceling the twin SU ( 2 ) anomaly requires an additional doublet neutral under 
the twin SU( 3): L , left-handed twin tau. Although not required for anomaly cancellation, 
introducing a twin right-handed tau t r allows the f to be rendered massive provided a 
Yukawa coupling 

C D y T BLf R . (10) 

As in the case of the twin bottom Yukawa, y T need not respect the Z 2 as long as y T <C y t . 
The twin neutrino may be rendered massive in the same way as the SM neutrinos; its mass 
plays no role in naturalness and is essentially a free parameter of the model. Finally, twin 
light-flavor (first and second generation) fermions are totally unnecessary for naturalness, 
as their Yukawa couplings too small to meaningfully disturb the Higgs potential, and are 
therefore absent in our minimal “Fraternal” model. 


C. Summary of Fraternal Twin Higgs model 


Thus we arrive at the ingredients of our minimal Twin Higgs model: 


1. An additional twin Higgs doublet and an approximately SU (4)-symmetric potential. 

2. Twin tops and a twin top Yukawa that is numerically very close to the SM top Yukawa. 

3. Twin weak bosons from the gauged SU( 2 ) with g 2 ( A) « g 2 ( A). 


4. Twin glue, a gauged SU(3) symmetry with g 3 (A) ~ < 73 (A). This gauge group is 
asymptotically free and it confines, at a scale A 3 that we will discuss in Section 


IV B 
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5. Twin bottoms and twin taus, whose masses are essentially free parameters so long as 
they remain much lighter than the twin top. 

6 . Twin neutrino from the twin tau doublet, which may have a Majorana mass, again a 
free parameter as long as it is sufficiently light. 

As an aside, we note that in contrast to a model with a perfect twin of the Standard 
Model, this model is cosmologically safe; with at worst one massless particle (the twin tau 
neutrino), and fewer degrees of freedom than the visible sector, the effective number of 
degrees of freedom during nucleosynthesis and recombination is very small. 

Not accidentally, the most crucial ingredients for the naturalness of the theory — a twin 
Higgs, twin tops, and twin glue — strongly resemble the ingredients of natural supersymme¬ 
try (Higgsinos, stops, and the gluino). But the key difference here is the likely existence of a 
new, confining gauge group in the minimal twin sector. Although twin glue does not impact 
the Higgs directly at one loop, its contributions at two loops make it a key component of a 
viable twin Higgs model. 


III. ELECTROWEAK BREAKING AND TUNING 


In this section, we study the effective potential of the Fraternal Twin Higgs model outlined 
above. We show how realistic electroweak symmetry breaking can be achieved, accompanied 


by a 125 GeV Higgs scalar, and estimate the tuning of couplings needed. In subsection III A 


we write down the effective potential for the full Higgs sector of the model at one-loop order. 
In subsection |III B[ we integrate out the heavier Higgs to get an effective potential for just 


the SM-likc Higgs. In subsection IIIC we determine that the degree of fine-tuning needed 
to have realistic electroweak breaking and SM-like Higgs mass is 2 v 2 /f 2 . Finally, in 


subsection III D we more fully justify the form of our starting effective potential in subsection 


III A, by showing that it is free from other types of fine-tuning. 
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A. Effective Potential 


The Twin Higgs effective potential is given, to good approximation, by 


Veff = A(|A | 2 + \B\ 2 — / 2 / 2) 2 
+ k(\A\ 4 + \B\ 4 ) 

+ p\A\ A + crf 2 \A\ 2 


167T 2 


[|^| 4 ln(Ay B ,V| 2 ) + |B| 4 ln(A 2 /y?|B| 2 )] , 


( 11 ) 

( 12 ) 

(13) 

(14) 


where for concrete estimates we will take the UV cutoff of the Twin Higgs theory to be 
A = 5 TeV. Lines 0 - ( |T3j ) represent the most general renormalizable tree-level terms 
consistent with the SM and twin gauge symmetries. Line ( fllj ) is just our starting point, 
Eq. (JTJ) , the subset of terms respecting the global £77(4) symmetry, under which (A, B) 


transform in the fundamental representation. Line (12) consists of the extra terms allowed 


by breaking £77(4) but preserving the discrete Z 2 global subgroup, A o B. (The other Z 2 
invariant, |A| 2 |£| 2 , is equivalent to this, modulo SU (4)-invariant terms.) Line (13) consists 


of the remaining extra terms which respect only the gauge symmetries. Line (14) is the 


dominant one-loop radiative correction that cannot be absorbed into a redefinition of the 
tree potential, due to a (twin) top loop. This is just Eq. ([3]) , where we have set y t = yp, 


in Section [TVj we will justify this as a good approximation because of the gauging of twin 
color. While the logarithmic cutoff dependence above can be removed by renormalization 
of k, we will keep the above form for later convenience in making estimates. 


B. Matching to SM Effective Field Theory 

Before fully justifying the above approximate structure for the effective potential, we will 
first work out its consequences, in particular matching it to a SM effective field theory at 
lower energies. This will allow us to choose the rough sizes of the different couplings needed 
for realism, and to then self-consistently check our approximation. 

Line 0 contains the dominant mass scale, set by /. We take the SU (4)-symmetric 
self-interaction to be modestly weak, A < 1. We seek to stabilize the vacuum such that 


/ = few x va, va = v = 246 GeV, 


(15) 
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so that the pseudo-Goldstone Higgs is primarily SM-like. We will also use this small hierarchy 
to work perturbatively in powers of v/f. 

We begin by studying the limit v/f = 0. The dominant potential, 0, then gives rise 
to a spontaneous symmetry breaking S77(4) —> SU( 3), where the breaking VEV is in the B 
direction, 

(B 0 ) = f/V2. (16) 


This breaks twin 577(2) gauge symmetry, with three of the seven Nambu-Goldstone bosons 
being eaten in the process. The remaining four real Nambu-Goldstone bosons form the 
complex A doublet, namely the SM Higgs doublet of electroweak gauge symmetry. The 
fourth uneaten B scalar is the radial mode of the potential, a physical heavy exotic Higgs 
particle with mass, 


m; 


•„ = V2A /. (17) 

Later, 0{y 2 j f 2 ) perturbations will mix these boson identifications to a small extent. 

The couplings p and a break the discrete Z 2 symmetry. The only other couplings in the 
theory that break Z 2 are the small SM hypercharge coupling and the very small Yukawa 
couplings to the light fermions (and their twin equivalents in the case of the r and b , also 
taken to be <C 1). Therefore it is technically natural to take p to be as small as rv_/ g\/(lQ-K 2 ), 
where gi is the SM hypercharge coupling. This is so small that we neglect it in what follows. 4 
We can however consistently take af 2 to represent an explicit soft breaking of Z 2 symmetry. 
Still, we will take cr < A so that line 0 dominates the mass scales as assumed above. 

At energies well below the heavy Higgs mass and heavy twin gauge boson masses, set by 


/, \A\ 2 + \B\ 2 is rigidly fixed at the bottom of the potential in line (11), 

f 2 

\B\ 2 = J —~ \A\ 2 . 


(18) 


Plugging this into the (A, B ) effective potential (and neglecting p as discussed above) gives 


4 The possibility of larger p offers alternate model building opportunities but we do not pursue them in this 
paper. 
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an effective potential for just the lighter A degrees of freedom, 


Veff(A) = 


2a — 2 k 




mi 


2k+ w ( ln < A2 M 2 w 2 ) + hi ( A2 /™ 2 ) - \ 

0(A e /f), 




( 19 ) 


where this equation is expressed in terms of the mass of the twin top, or “top partner”, 

/ 


= y t 


V2 1 


( 20 ) 


in analogy with the top quark mass. 


Eq. (19) has the form of a SM effective potential, with the tree-like \A\ 2 , | A 1 terms, as well 
as a top-loop induced \A\ A ln(A 2 /?//|A| 2 ) term. Successful electroweak symmetry breaking, 
(A 0 ) = 246GeV/-\/2, and a physical Higgs mass of 125 GeV, can therefore be arranged by 
tuning the \A\ 4 coefficient using k, and the \A\ 2 coefficient using a. To estimate the tuning 
involved, we can neglect the modest In |A| modulation of |A| 4 , and just set the logarithm to 
its expectation value, 




eff 


2(T ~ ^ ~ S ( hl(A/mf:) \ 

Asmv 2 \A\ 2 + X sm \A \ 4 , 


9 

Vfli 

IF 1 ' 41 ' 


2k + wl ln ( A M~) + - - A 




( 21 ) 


where rn t = yt(Ao). For realistic electroweak scale and physical Higgs mass we require 
Asm ~ v = 246 GeV. 

Noting that In (rrif/mt) — | = ln(few) — | ~ 0(1), we will neglect this term relative to 
the ln(A/rrq)-erihanced top-loop contribution to the Higgs quartic, so that 

3w 4 

A sm ~ 2k + ^ ln(A/mf). (22) 

With t/ f ~ 1, the top-loop induced quartic coupling is already close to the required value, 
A sm ~ 1/8. A higher-order and renormalization-group improved analysis, as in studies of 
analogous corrections in supersymmetry, is expected to reduce this radiative correction, the 
central feature of which can be captured by using a top-Yukawa coupling renormalized at 
several TeV, where yf ~ 1/2. In any case, the rough sizes of the radiative corrections in 
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the Twin Higgs theory are comparable to the realistic value of A sm, and so with a n of 
comparable magnitude we are able to successfully and naturally fit the observed physical 
Higgs mass. 

As an aside, it is instructive to compare the form of the logarithmically divergent (twin) 


top-loop contributions in Eq. (19) with the analogous contributions from top/stop loops in 
weak scale supersymmetry (for large tan /3 and small stop mixing). As in supersymmetry, the 
quadratic divergence in the top-loop contribution to the Higgs potential has been canceled by 
a top partner mechanism. However, there remains a logarithmically divergent contribution 
to the Higgs mass-squared. This has precisely the same form as in supersymmetry with 
the replacement, m s t op —>• rnp We see there are also logarithmically divergent contributions 
to the Higgs quartic self-coupling from the (twin) top loop. Again, there are analogous 
contributions in supersymmetry of the same magnitude, but the stop contribution has the 
opposite sign due to its Bose statistics. Thus in supersymmetry the logarithmic divergence 
cancels out here, but there is still a finite logarithmic enhancement factor of In {m 2 ^ /rn\) to 
the Higgs quartic correction. Unlike the MSSM where the tree-level Higgs quartic is set by 
electroweak gauge couplings, here we have an unconstrained tree-level quartic contribution, 
2k. (However, k also appears in the Higgs mass term, so it will affect electroweak tuning, 
as discussed below.) 


C. Estimating Electroweak Tuning 


Electroweak tuning involves the quadratic terms of the potential. Fortunately, the same 
combination of parameters that dominates the quartic self-coupling, in the second line of 


Eq. (21), also appears in the quadratic terms of the potential, in the first line of Eq. (21). 


We will drop the 1/4 term relative to the In A-enhanced part of the top radiative correction. 

Then, matching the quadratic terms to the SM form on the last line, we have 

2 

777 /- 

A smV 2 ~ ( A sm ~ 2cr)—A . (23) 

yf 


The degree of electroweak tuning of a needed to achieve this is therefore 

A smv 2 


Electroweak — Tuning 


A SM(m 2 { /y?) 

2 m 2 t _ 2v 2 

= 7 ^' 


m- 


(24) 
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For example, for / ~ 3v (rrif ~ 500 — 600 GeV), this corresponds to a very mild 20 percent 
clectroweak tuning. 


D. Twin Higgs Effective Potential Approximation 


Finally, we justify our starting approximation for the Twin Higgs effective potential upon 
which our analysis of electroweak breaking and tuning is based. To do this, we examine each 


term of Eqs. (11 - 14) and ask whether its coefficient is radiatively stable, and also whether 


this coupling itself significantly radiatively corrects other couplings. 

The tuning needed to have / smaller than the cutoff A is determined by the leading 
radiative corrections to the quadratic terms in line (JTTj) . We have taken the SU (4)-symmetric 
self-interaction to be weakly coupled enough that the leading quadratic radiative corrections 
come from the (twin) top loop and the large Yukawa coupling, 


V {1 - loop ) D —-|/, 2 A 2 (\A\ 2 + \B\ 2 ) . 

8tt 2 v ' 


(25) 


The other couplings that can similarly contribute at one loop are the electroweak gauge 
couplings, as well as k and p. The electroweak couplings are subdominant to the top 
Yukawa coupling as reviewed earlier. We have estimated that k 0(X SM ) (Eq. (|22j) and 
ensuing discussion), and loops using the SM Higgs quartic coupling are also considerably 
subdominant to top loops. Finally, we have consistently chosen p to be very small as 


explained earlier. Therefore, Eq. (25) dominates radiative corrections to A/ 2 and gives a 
tuning to keep / A of 

A/ 2 8 tt 2 A/ 2 


/ — tuning 


(26) 


sfi» 2 A 2 3yfA 2 

If we take for example / ~ 3v ~ 750 GeV, we see that there is essentially no tuning. 

We are taking k to be comparable to the leading radiative corrections from the top 


loop (Eq.(22) and ensuing discussion), so this Z 2 -preserving but SU (4}-violating coupling 


is radiatively stable and natural. As explained earlier, we can naturally take the hard 
breaking of Z 2 given by p to be negligibly small, erf 2 represents a soft breaking of Z 2 , with 
cr ~ Asm/2 (Eq. (23) and ensuing discussion), and is clearly natural without other sources 


of soft breaking. To tune to realistic electroweak symmetry breaking we saw that we needed 
a ~ Asm/2 ~ 1/16, therefore we can easily have cr/ 2 <C A/ 2 to ensure that line ( |IT| ) indeed 
dominates the mass scales of the Twin Higgs potential as our analysis presumes. 
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With all couplings in our theory being perturbative, our one-loop analysis suffices for 
demonstrating the successful matching to a realistic SM effective field theory and for esti¬ 
mating the tuning required. 

IV. FRATERNAL COLOR 

We now discuss the dynamics of twin color in more detail. We first calculate how close the 
twin and visible sector color couplings must be to preserve naturalness. Then we estimate the 
confinement scale of twin color and discuss the associated twin hadrons, including glueballs 
and quarkonia. 


A. Perturbative Considerations 

As we have seen above, the twin Higgs mechanism for naturalness, at one-loop order, 
requires the top and twin-top Yukawa couplings to be very nearly identical close to the 
cutoff. Here we show that when QCD effects are taken into account at two-loop order, 
naturalness favors having a twin QCD, with a gauge coupling similar to that of QCD near 
the cutoff A. 

To see this, consider the one-loop RG analysis for the dimensionless Wilsonian running 

2 / \ 

mass-squared parameter of the SM-like Higgs, x(g) = "y : 

dx = o r | 3 (^ 2 - y't ) . 
d In g ' 2tt 2 

dyt _ 9 yf y t g\ dy t _ 9 $ y t g\ 

d\wg 327T 2 27r 2 ’ d In g 327T 2 27 t 2 ’ 

dgs = 7g| dg 3 = 29^| 

d In g 167T 2 ’ d In g 487r 2 

where g 3 is the SM QCD coupling and g 3 is the twin QCD coupling. Note that g 3 and g 3 
run differently because the SM has six QCD quark flavors while the minimal model has just 
two twin QCD flavors. If we neglect the running of yt, yt, ^ 3 , § 3 , then the solution to the first 
of these equations is simply given by 

+ %, 4~/‘ 2 V - A 2 ) • (28) 

In running down to the weak scale or physical Higgs mass, /1 C A, we see we have merely 
matched onto the one-loop result of Eq. (|5]) , which stresses the importance of having y t and 
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y t be very nearly the same. But taking account of the running of yt,yt, fiq, <?3 by solving all 


of Eqs. (27) then gives an RG-improved result, which allows us to explore the role of twin 
color in maintaining y t ~ y t as the couplings run. 

The simplest calculation arises by seeing what happens when we have no twin color, 
setting (73 = 0. The danger to naturalness is then that the SM QCD coupling will give 
yt(l-i) a different evolution to ytfg) even if they happen to coincide at /i = A, and that this 
deviation will feed into m 2 . To focus on just this effect we will drop the (^-independent 
terms in the y t , yt /3-functions, since these have identical forms for the twin and SM sectors. 
(Keeping these effects would be subleading in the running of m 2 compared to those of the 
(^-dependent term if we start with y t ( A) ra y t ( A).) As a final simplifying approximation we 
will drop the running of g 3 itself, because its /3-function is small near A and because we will 
see that the running in m\ is dominated near A. With these approximations, and working 


to first order in y t — y t , the solution to Eqs. (27) is given by 

^yhl , 3 (ti-y\ - 


m h (/ 1 ) « m h { A) + 


87T 4 


+ 


47T 2 


2\ . 21 / /*\ 

h ) + ln W a ) • 


(A- 


(29) 


Running down to low scales, /i<A, 

™h,m ~ " 4 ( A ) + 


3Vt9l , 3 {y 2 t 
87r 4 


+ 


vt. 


47T 2 


A 2 


(30) 


In the above equations all the dimensionless couplings are evaluated at A. 

We see that even if y t ( A) = y t ( A), just the running from SM QCD has led to a quadratic 
divergence which would require a fine tuning of the counter-term m\{ A) to get the physical 
Higgs mass, 

0 

TTl 

Two loop fine tuning ~ - ~ 0.25. (31) 

8tt 4 Vi 

Of course, it is not reasonable for two couplings, y t , yt , that run differently to be exactly the 
same at A. This would be a fine tuning in itself. A technically natural estimate is that they 
differ by at least the running of y t due to QCD over an e-folding of running, 

MA)-S«(A)|>dA»(A). (32) 


In other words, we expect a comparable quadratic divergence in m 2 IR from the splitting 
in Yukawa couplings and from the explicit Oig^yf) divergence. Thus a better estimate of 
fine-tuning in the absence of twin QCD is < 10 percent. 
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The tuning due to QCD at two-loop order will clearly become very mild if we do include 
twin QCD with g 3 (A) ~ ( 73 (A). The two couplings will run differently in the IR due to 
the different number of SM and twin generations, but quadratic sensitivity to A will be 


determined by the UV couplings. Roughly, Eq. (30) is then replaced by 

-3y?(AM(A)-9l(A)) , 3(y t 2 (A)-#(A)) 


™h,IR ~ "%( A ) + 


+ 


A 2 


(33) 


87 t 4 4n 2 

Therefore, if g 3 (A) and g 3 (A) agree to within even 15%, enforced by an approximate Z 2 , 
the tuning in m 2 will be a mild ~ 30%, comparable to that in Eq. (24). This estimate 


combines the 0{g\yf) contribution from running and the threshold correction in Eq. (32) 
in quadrature to reflect the unknown relative sign of the threshold correction. We will take 
this to be the case in what follows and study the effects of twin confinement from this twin 
QCD sector on the twin spectrum and phenomenology. 

Finally, of course even g 3 and g 3 run differently because the particle content of the minimal 
Twin sector differs from that of the SM, so one may wonder how close they can naturally 


be at A. However, the analogous estimate to Eq. (32) is 

|ga(A) -&(A)| ^ g§(A) 
g 3 (A) ~ 6 vr 2 

which is easily consistent with the requirement of naturalness discussed above. 


(34) 


B. Fraternal Confinement 

We now determine the confinement scale A 3 of the twin SU (3) gauge interaction, as this 
governs the infrared phenomenology of the twin sector. 

If all fermions carrying twin SU (3) quantum numbers are much heavier than the confine¬ 
ment scale, the infrared physics is that of pure SU( 3) gauge fields, and the lightest states in 
the confined twin sector will be glueballs, whose rich spectrum includes states with different 
angular momentum J, charge conjugation C and parity P. As shown in lattice studies of 
SU (3) pure glue [441145] , at least a dozen glueballs are stable against decay to other glueballs, 
and the lightest, which we will refer to as G 0 +, Las J PC — 0 ++ . Lattice data provides the 
ratios of these glueballs’ masses to each other and to the confinement scale. To determine 
the physics of the twin sector thus only requires us to compute A 3 . 

Once the twin b becomes sufficiently light that twin glueballs and twin bottomonium 
states (which we will refer to generically as “[56]”) have comparable masses, the situation 
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FIG. 3: The confinement scale A 3 of the twin SU{3) coupling given fractional variations in 53 
and y b at the cutoff A for the minimal Twin Higgs (dependence on y b is negligible). Here we take 
A = lOrht and / = 3v. The mild kinks are due to the b threshold. 


becomes more complex. We will not explore this regime carefully in this paper, leaving its 
details for future study. However, the calculation of A 3 and of glueball masses given below 
still applies approximately. 

The twin and SM SU(3) couplings are similar at the cutoff A, but the twin sector has 
fewer quark flavors, faster running (i.e. a more negative beta function), and therefore a 
modestly higher confinement scale. This is illustrated in Fig. [3j which shows the strong 
coupling scale A 3 of twin QCD as a function of the variation dg 3 between SM and twin QCD 
couplings at the cutoff as well as the value of y b relative to y b . Note that for g 3 ~ g 3 , A 3 is 
typically one to two orders of magnitude above that of QCD, with weak dependence on y b 
through its impact on the twin QCD beta function . 5 

We may now estimate the mass scale of twin glueballs. Using lattice estimates of the 


At two loops, we define the MS confinement scale A 3 via 


Af s 

y 


= exp - 


2& 0 ff!(/d 


(bogl(y)) 


-b 1 /2b 2 0 


1 + T-gl(y) 

° 0 


bi/2b 2 0 


(35) 


where AbU are the one-loop and two-loop twin QCD beta functions respectively, and <73 (y) is understood 
to be the MS coupling. 
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glueball mass spectrum in units of the inverse force radius [44] 05J and the zero-flavor SU (3) 
MS confinement scale in units of the inverse force radius [46], we find the mass m 0 of the 
G 0 + glueball to be related to the strong coupling scale via m 0 = 6.8A^ IS . The physical scale 
mo determined by running couplings down from the cutoff carries a combined uncertainty 
of 0(10%) from the lattice estimates, primarily due to uncertainty in the inverse force 
radius ffl. Given the mass of the Go+, the masses of higher glueball excitations in terms 
of mo are known to good precision. The next highest states in the glueball spectrum are 
well separated from the G 0 +, with the closest states being the G 2 +, a 2 ++ glueball with 
m 2++ ~ 1.4 m 0 , and the Go-) a 0 _+ state with m 0 _+ ~ 1.5 m 0 [44J 45]. It appears also that 
there is a second stable 0 ++ glueball, the G' 0+ , with m ~ 1.8 mo m- The twin glueball 
spectrum is illustrated on the left-hand side of Fig. [4j 

Meanwhile, the twin bottomonium states form a rich spectrum, whose lowest lying states 
are narrow if rapid decays via twin glueballs are inaccessible. As is familiar from SM quarko- 
niurn data, and as sketched on the right-hand side of Fig. [4| the spectrum includes towers 
of 0 +_ , 1 and j ++ (j = 0,1, 2) states, which by analogy we call fj, T and x.y The lowest- 
lying twin quarkonia states have masses of order 2mg, with mass splittings -C 2 m b (as long 
as m~ b A 3 .) However, since there are no other light twin quarks, there are no “open twin 
bottom” mesons analogous to the SM’s B u , Bd mesons. Thus the towers of narrow quarko- 
niurn states extend much further up than in the SM, potentially up to a scale of order 2mo, 
as sketched in Fig. [4j or m 0 + 2mg. 

The reader should note that Fig. [4] is only illustrative, and must be interpreted with 
caution. Its details change dramatically as one raises or lowers the quarkonium masses 
relative to the glueball masses, and it omits the many narrow higher-spin quarkonium states, 
along with various other phenomenological details. 


V. TWIN HADRON PHENOMENOLOGY 

Thus far we have seen that viable Fraternal Twin Higgs models include twin glue, with 
couplings that favor confinement roughly an order of magnitude or so larger than the SM 
QCD scale, producing relatively light twin glueballs and/or twin quarkonia. Both twin 
gluons and twin quarks are connected to the Standard Model via low-dimensional portals, 
and this can lead to observable and even spectacular twin hadron phenomenology. As in 
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FIG. 4: Sketch of the twin hadron spectrum in the regime where mo < 2mr < 2mo- In addition to 
the Go+, of mass mo, about a dozen other glueballs, with mass splittings of order mo, are stable 
against twin strong decays. Numerous twin bottomonium states, including a tower of 0 ++ states 
X, are stable against twin strong decays. The circled Go+ and G' 0+ glueballs, and potentially the x 
quarkonia, can dominantly decay via annihilation through an s-channel off-shell Higgs to the SM. 


Folded Supersymmetry {4j, where twin glueballs also arise, we thus find a connection between 
dark naturalness and twin hadrons. In our case, this connection manifests itself as new and 
exciting opportunities for discovery at the LHC. 

The model’s phenomenology changes significantly as we move around in the parameter 
space, and in most regions it is rather complicated. But the most promising and dramatic 
LHC signals arise even in the conceptually simplest region, namely where mg > \rrih (i.e., in 
Fig. [5] below, the part of region A above the dashed line). In this case the main phenomenon 
is that described in Fig. [lj with twin gluons produced in h decays and hadronizing into 


twin glueballs, including the Gq+. The Gq+ lifetime is discussed in Section V C Eqs. (39) - 


(41); the (perturbative) production rate for glueballs is discussed in Section VE, Eq. (45), 


with nonperturbative subtleties described in Appendix B 1 The reader seeking to avoid 
becoming lost in details at a first reading may wish to focus merely on this simple scenario, 
in which case Section V D, the later portions of Section V E , and Appendices [A] and B 2 
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Mass of Go- glueball (GeV) 


FIG. 5: The parameter space of the model in terms of the masses of the lightest glueball Go+ 
and the lightest quarkonium fj. In region A, only glueballs are produced; in region B, the relevant 
quarkonia decay to glueballs; in region C, glueballs are either not produced or decay to quarkonia, 
so only quarkonia appear in the final state; and in region D there are both metastable glueballs and 
metastable quarkonia, with the potential for mixing. Solid lines indicate kinematic boundaries. 

may be omitted. 


A. Kinematic Regions 

Before we begin, it is useful to parameterize the theory through mo and m,- ? (as well as /) 
in place of < 73 , y^. Here fj is the lightest [bb] state, lying slightly below the lightest x state. We 
can then divide the parameter space of the model into four qualitatively different kinematic 
regions, shown in Fig. [5j 

• Region A: > 2mo, iTth < 2and rrih < m 0 + rrif n so that h can decay to twin 

glueballs but not to twin bottomonium. 

• Region B: nih > mo + m^ and rrifj > 2mo; here h can produce twin bottomonium, but 
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all x states (and most or all other [bb] states) decay eventually to twin glueballs. 

• Region C: Either mo + rn^ > mj l > 2m^, in which case only twin bottomonium can be 
produced in h decays, or m 0 > 2m,), so that any produced glueballs decay rapidly to 
twin bottomonium. 6 Only bottomonia appear in the final states. 

• Region D: mh > mo + m^, < 2mo, mo > 2m^; here both the lighter twin bot¬ 

tomonium and glueball states are metastable, and either or both may appear in final 
states. Mixing between the two classes of states can be important in this region. 

As we will see, visible decays are typical, and displaced decays are possible, in regions A, 
B and D. They may or may not be absent (or very rare) in most or all of region C; this is 
model-dependent. 


B. Couplings to the Visible Sector 


In this model, the two portals of greatest importance will involve the dimension-five 
operator Abb and the dimension-six operator A A G I1V G I1U , where recall A is the SM-like 
Higgs doublet. After SM electroweak symmetry breaking and twin confinement, the first 
operator causes mixing between the SM-like Higgs and twin x quarkonia, while the second 
operator causes mixing between the SM-like Higgs and twin Go+ glueballs. The effective 
dimension-5 coupling A'Abb originates from the twin bottom Yukawa ybBQd ; applying (18) 
gives the leading interaction 

£5 = ~yb^i=r;bb, (36) 


72/ 


which after electroweak symmetry breaking yields C D —^'jhbb. This is just the v/f- 
suppressed coupling of the SM-like Higgs to twin fermions, which gives rise to mixing between 
h and x quarkonia. 7 


6 Not all of region C is physically meaningful. As mo rises, so does A 3 , and mfj always remains heavier than 
the confinement scale. 

7 No other portals are relevant for h decays. Below the heavy h scale, and in the absence of a twin U( 1), 
there are no dinrension-4 portals; however we do discuss the effect of these portals in Section VI C| and 
Appendix [E] A dimension-5 portal involving SM and twin neutrino mixing is not induced or required in 
our model and may easily be too small to have any measurable effect. 
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While the value of the h — \ mixing can vary depending on the (unknown) value of the 
twin bottom Yukawa, the h — G 0+ coupling is necessarily generated by the basic ingredients 
of the minimal Twin Higgs. In a manner entirely analogous to the Standard Model hgg 
coupling, loops of twin tops generate an effective coupling between the twin Higgs doublet 
B and twin gluons; after 577(4) breaking and SM electroweak symmetry breaking, this leads 
to a coupling between twin glue and the SM-like Higgs. The effective coupling between the 
twin Higgs doublet B and twin glue takes the usual form 


a, - ~ f'B^B\ 

jC 6 = In ' 


12tt 


(iv ^ a 


p J 


(37) 


and applying (18) generates the corresponding coupling to A* A, which after electroweak 

(38) 


symmetry breaking leads to 

Here the couplings are largely fixed by naturalness considerations. 


r ^ r^ 

6tt//^“ ' 


C. Glueball Decay 


Once produced, twin glueballs can decay into kinematically available final states. As¬ 
suming there are no light quarks in the twin sector, the only potentially available decays 
are into light Standard Model fermions via the (off-shell) SM-like Higgs h or into the twin 
lepton sector via the heavier Higgs h or twin Z. The decay Go+ —>• h* —>• YY, where Y are 
light SM fields, provides a visible signal. This process was studied in [23J in the context of 
a similar Hidden Valley model [21] . with a similar Higgs portal for decays of the glueballs, 
and with a production portal induced by “quirks” [13 122]- Tlie width for G 0+ YY is 

rG 0+ ^yy = ( -2^ r f^yy( m o) • ( 39 ) 


V 6v Tf 2 (m 2 h -m 2 0 )J h ^ YYK 

Here T^ yy (mg) is the width of a SM-like Higgs of mass rrih = mo and fo is the G 0 + decay 
constant; from the lattice we have 47rb3fo = 3.06mQ [[45] . This is the dominant decay mode 
of twin glueballs in the minimal model; decays into twin sector leptons are subleading, due 
to suppression factors of (mh/mf) A (f /v) 2 for decays via h and by an extra (u//) 2 mixing 
factor if via h. 

Provided that the Go+ decays primarily into SM final states, we can determine its lifetime 
tq in terms of mo and /. This exercise is particularly straightforward for glueballs much 
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cr(0 ++ )[log 10 (m)] 



FIG. 6 : Decay length ctq of the Gq+ state in log 10 (meters) as a function of mo and /. 


lighter than the massive Standard Model gauge bosons, for which the factor r^fyy(mo) 
scales linearly with glueball mass because decays into longitudinally-coupled modes are 
doubly off-shell. Then in this regime we have simply 


r - i.i x nr 17 GeV x 


m 0 

10 GeV 


7 / 750 GeV 

\—J— 


4 


? 


(40) 


valid for 2 m T < m 0 < mw■ This corresponds to a decay length ctq of approximately 


CTo ~ 18 m x 


10 GeV\ 7 / / 

m 0 ) \750 GeV 


4 


(41) 


Note mo ~ 10 — 50 GeV is a central range of values for the glueball mass, given (a) the 
relative factor of ~ 7 relating A 3 to mo and (b) the higher confinement scale of twin QCD 
given a reduced number of twin sector fermions, see Fig. [3j This is a tantalizing result from 
an experimental perspective; it implies that the twin sector glueballs give rise to displaced 
decays on the length scale of the LHC detectors. Even for average decay lengths greater than 
the scale of the detectors, the large number of h bosons produced at the LHC means that 
an appreciable number of distinctive displaced decays can still occur within the detector. 
These decays involve a mix of final states with relative rates corresponding to the decay of 
an SM-like Higgs of mass m^ = m 0 — primarily bottom quarks, if kinematically accessible, 
as well as tau pairs and gluons. 
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We can improve upon the above parametrics by computing the SM Higgs width as a 
function of mass using HDECAY [0Sj. A precise plot of the decay length cr 0 as a function 
of m 0 and / is then shown in Fig. [6] As we can see, the decay length consistent with 
naturalness in the minimal model is typically in the range of millimeters to meters, with a 
jump toward the kilometer scale for mo < 2 mb- 

Note, however, that given the strong dependence of the glueball decay length on the 
glueball mass, cto oc tHq, even modest uncertainties in the value of mo may substantially 
influence this estimate. For example, the 0(10%) uncertainty on the value of mc 0+ inferred 
from lattice data translates into a factor of two variation in the value of ctq. 

The second 0 ++ glueball G' 0+ decays similarly to the first, with a shorter lifetime, and 
may have prompt or displaced decays. Others such as the 2 ++ may decay via an off-shell 
Higgs, e.g. G 2 + —> Go + h*. At small masses their lifetimes are very long |23j, and they will 
appear only as missing energy in h decays. At larger m 0 they cannot appear in h decays, 
though they might possibly appear as displaced vertices in decays of the heavy twin particles 


Z or h. 


The above calculations apply in regions A and B but are irrelevant in region C, where in¬ 


stead G*o+ —> [bb] [bb]. They are only part of the story in region D due to glueball-quarkonium 


mixing, which we briefly address at the end of Section V D 


D. Bottomonium Decay 


The twin bottomonium states shown in Fig. [4] resemble the SM’s towers of charmo- 


nium/bottomonium states, yet have unfamiliar decay modes. We relegate most details to 
Appendix [Aj but summarize key points here. Like the G 0+ , H ie \ states can potentially 
decay via annihilation to SM states through an s-channel Higgs. However, if twin neutrinos 


(and/or taus) are light, they can also decay to spin-one T states by twin weak interactions, 
X —>■ T££, where £ — is, r. If allowed, this decay competes with the decay via the Higgs, and 
the competition is very sensitive to modeling of the [bb] states. 


We can use non-relativistic quantum mechanics to give a rough leading-order estimate of 
the width of the y state to the SM: 



( 42 ) 
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where R{r ) is the state’s radial wave function. Here R'( 0) appears because x is a p-wave 
state. As justified in Appendix [A] we take the approximation of a linear confining potential 
(with slope a ~ 4Ag, following calculations of [441 49j) as a starting point, ignoring the 
Coulomb potential and important relativistic corrections. Our estimate is 


x->yy 


2 x 10 


-3 


11/3 10/3 


v\ m x m 0 
f) v 2 m h (m 2 h - ml) 2 h ^ YY 


(m h ) 


(43) 


for the lowest x state. 

Meanwhile the twin weak decay to Til proceeds via an off-shell Z, through a dipole 
transition (analogous to Xb —^ T 7 in the SM). Assuming twin neutrinos are massless we 
roughly estimate 


a 2 \ m x ~ m r J 


x^ruo 4?r ( ms/ g)2 m | 


(44) 


where f3 is the typical velocity of the b in this state. The extreme dependence of this width 
on the x — T mass difference makes any estimate of lifetime and branching fraction highly 
uncertain. 

Using these formulas as a guide, however, we can qualitatively summarize x phenomenol¬ 
ogy in regions C and D: 

• If twin neutrinos and taus are not light, then one or more low-lying x states may decay 
promptly at high mass and displaced at low mass. An approximate formula for the 


lifetime is given in Eq. (A3). 


If twin neutrinos and/or taus are light, then 


— I 11 region C, twin weak decays dominate, making all x decays invisible. 

— I 11 region D, the lowest-lying x state may decay visibly with a substantial branch¬ 
ing fraction, especially for m x > ?n 0 . 

— Also in D, the lowest-lying x decays may be displaced for m x > mo and m x < 
30 GeV or so. 


Recall also that in region D there can be mixing of 0 ++ glueballs and quarkonia. This effect 
is probably of greatest important for m 0 > m x , where G$ + is nested within the x tower. For 
mo ~ the quarkonia widths tend to be larger, so the mixed states tend to inherit their 
properties, with lifetimes somewhat shorter than in Fig. [6] and the potential for significant 
invisible decay fractions. 
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E. Twin Hadron Production 


We now discuss twin hadron production via h decays. 8 Twin hadronization is complicated 
and quite different from SM QCD, so it is not possible to make reliable estimates as to what 
fractions of twin glueballs and bottomonia produced are the golden 0 ++ states that we can 
hope to observe. But the overall rate for h —> twin hadrons, which occurs through h —> gg 
(if m 0 < mh/2) and through h —> bb (if m b < m^/2), can be estimated, at least roughly. 

We begin with a perturbative analysis of h decays to the twin sector; non-perturbative 
effects, which can qualitatively change the story, will be discussed later. We will refer to the 
h —>■ gg production of glueballs (as shown in Section [TJ Fig. [l]) as the “irreducible process”, 


since, proceeding via the interaction (38), it is independent of m b to a first approximation. 


It is the only relevant process for m b > rrih/ 2, and also dominates for yb < 0.13yb- The 


partial width for the irreducible process follows from (38), 

- 9 \ 2 

«3 v 


r(fc -> gg) ^ I — -p ] T(h-tgg), 


(45) 


leading to a perturbative expectation for the branching ratio of the Higgs to twin glueballs 
of the order 0.1% for / = 3v. 

Similarly the decay h —$■ bb can generate [bb] states in regions B, C and D, and also 
glueball states which may be produced immediately along with the [bb] ( e.g. h —* xG 0+ ), or, 
if kinematically allowed, in radiative decays ([bb] —>• [bb]' + glueball) or via bb annihilation 
( [bb] —> glueballs). In region A, [bb] states are inaccessible, but still this process can give 
non-perturbative enhancement to twin glueball production. 

At / ~ 3v the perturbative h ^ bb width equals the irreducible width for yj, ~ 0.13 yb and 
grows as y b , enhancing the irreducible rate by about ~ ^(yb/Ub ) 2 ■ This branching fraction 
becomes so large that the perturbative h —>■ bb rate would be inconsistent with current 
Higgs measurements for y b > 1.25 yb (at least for / ~ 3v] larger values of y b are allowed as 
/ increases, as we will discuss further in Section VID) — unless of course m b > rrih/Z, in 
which case only the irreducible process remains. 

In sum, within perturbation theory, 

• The high m b region — the upper portion of region A, where only glueballs are kinemat- 


Twin hadron production via the heavy h is discussed in Section VIC 
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ically allowed — manifests only the irreducible process, with Br {h — y twin glueballs) 
~ 0.1% for / ~ 3v, and decreasing like f~ A . 

• The low m b region, including low-mass portions of regions B, C and D, has a substantial 
h -X bb rate, with branching fraction increasing to ~ 10% near m b ~ 15 GeV (for 
/ ~ 3u). 

• Between these two regions the model is ruled out if we rely on perturbation theory. 
However, as we will immediately see, these perturbative considerations are not always 
applicable. 

As in SM QCD, where e + e _ —> hadrons is strongly modulated by resonances and other 
non-perturbative effects below a couple of GeV, non-perturbative effects are potentially 
substantial in h decays when kinematics restricts the number and types of hadrons that 
can be produced thereby. We discuss this in some detail in Appendix [Bj but here we just 
summarize our most important observations. 

• At moderate to large m 0 , higher 0 ++ glueball resonances (and gaps between them) can 
enhance (and suppress) the irreducible rate. A very conservative estimate is that, in 
the relevant kinematic regime, the suppression can be no worse than 1/10, and is likely 
less severe. Thus the irreducible rate is at least 10 -4 , and usually larger, throughout 
the parameter space. 

• At moderate to large m b) similar enhancement and reduction may occur from excited 
X states. If mh lies between two narrow x resonances, then there can be substantial 
reduction relative to the perturbative h -x bb prediction. As we discuss in a preliminary 
fashion in Appendix [Bj it appears that wherever we can calculate y widths and mass- 
splittings with order-one confidence — relatively large m b and relatively small mo 
the suppression between resonances is always enough to invalidate the perturbative 
exclusion. At larger mo (for fixed rrijf) the widths tend to be even narrower, the 
splittings larger, and the suppression somewhat greater, though not calculable with 
current methods. Due to this non-perturbative suppression, we are confident this 
model is not universally excluded when m b lies near and somewhat below m^/2. In 
fact we believe, for various reasons, that the model is still allowed in quite far into the 
lower m b portion of region A, and probably also survives in an upper portion (as well 
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as the lower portion) of region D, but much more work is needed to put this suspicion 
on a solid footing. 

In summary, we expect the branching fraction for h — y twin hadrons to be at least 10 -4 
and, in most of parameter space, much larger. Final states may include various combinations 
of twin glueballs and bottomonium. Although the rates for production of individual types 
of hadrons cannot be calculated, we expect a substantial fraction of the twin glueballs (and 
perhaps bottomonia) to be Go+ and \ states that can decay, possibly displaced, by mixing 
with the Higgs. While 1.25 m^ffv) < m~ b < ro/,/2 is excluded within perturbation theory by 
existing measurements of the Higgs, this part of parameter space is likely a patchwork of ex¬ 
cluded and allowed regions, whose details we cannot calculate without full non-perturbative 
information. In the allowed regions, the twin hadron production rate is always substantially 
enhanced beyond the irreducible rate. 

A qualitative summary of the phenomenology in the various regions is given in Fig. [7} 
We have conservatively assumed that twin weak decays of bottomonium are kinematically 
allowed, making decays in region C invisible, and making \ visible only in some limited 
portion of region D; if instead the twin weak decays are forbidden, \ decays are visible 
throughout C and D and displaced at low mass (see Appendix [Aj) In the next section we 
explore these and other signals more carefully. 


VI. LHC PHENOMENOLOGY 


By far the most spectacular signal that can arise from our minimal Twin Higgs model 
is the displaced decays of twin glueballs and quarkonia. We describe the phenomenology of 
this signal, as it arises from h decays, in Section |VI A[ further details on search methods 
are given in Appendix 0 If no displaced decays are observable, h decay signals may be 


challenging for LHC but certainly accessible at a future e + e“ collider (Section VIB). A 
brief discussion of how a twin hypercharge U( 1) would affect the phenomenology is given 
in Appendix [e} Section 


VIC covers signals from a heavier Higgs h. Finally, Section 


VID 


explores the order-(r;//) 2 effects on SM Higgs production rates. 
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FIG. 7: A qualitative overview of the phenomenology, for / = 3v, in the various regions of param¬ 
eter space; see Fig. [5} Details are explained in subsequent sections. Solid lines indicate kinematic 
boundaries. Common final states are indicated in italics. At low glueball mass, decays of the Go+ 
are displaced; see Fig. [6] Here it is assumed that there are light twin leptons, so one x state is 
visible, and even displaced, only in small regions; otherwise x decays visibly throughout regions C 
and D, and is displaced at low mass. 

A. New Higgs Decays With Displaced Vertices 

The branching fraction Br {h —» twin hadrons) > 10 -4 everywhere that it is not kinemat¬ 
ically forbidden. Because the number of Higgs bosons produced at LHC in Run II will be 
of order 10 7 , and because displaced vertices are spectacular signals when identified, these 
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numbers represent a very promising opportunity. Already hundreds or thousands of events 
with displaced vertices may have been produced, though in many parts of parameter space 
they would clearly have evaded existing LHC Run 1 searches [22H2HI EE EJ]- 

As Fig. [/[suggests, the model exhibits a great diversity of displaced vertex phenomenology. 
Rather than address the full story here, we mainly discuss the regions with the simplest 
phenomenology. These regions, it turns out, produce most of the possible Higgs decay 
signatures, and are therefore sufficient to motivate the most important searches, which are 
sensitive to effects in more complicated regions that we will not discuss in detail. 

The simplest region is the portion of region A with mg > m^/2, where the irreducible 
rate applies and only glueballs can be produced. As we move across this region from large 
to small mo, taking / ~ 3v, we find the following twin hadron phenomena: 

• For m 0 > 40 GeV, h —> G 0+ G 0+ dominates and G 0 + decays are prompt. 

• For 10 GeV < m 0 < 40 GeV, the G*o+ decays are displaced; decays to other glueballs, 
and to higher multiplicities of glueballs, become more common for smaller m q. The 
decay h —>■ G 0+ G" 0+ is of particular note, since the G' {)+ decays visibly. 

• Below about 10 GeV the Go+ lifetime is so large that decays in the detector are rare. 
This is partly compensated by higher glueball multiplicity per event. 


Consider next small mg < 15 GeV. In region B, where mo is small, the irreducible process 
produces the same phenomena as in region A, but h —> bb enhances the rate for twin hadron 
production, leading (via bottomonium decay to glueballs) to > 3 twin glueballs per event. 
In the low-mg portion of region C, the glueballs instead decay to bottomonium, leading 


to states that may all be invisible; alternatively (see Section VD) final states may include 
prompt or displaced \ decays. In region D, where mixing may be important, and where 
X —> SM has a larger branching fraction, an even richer set of final states is possible. 

With all of these different subregions with different phenomenological details, many of 
which cannot be calculated, one may rightly worry that experimental coverage of this model, 
and others like it, will be extremely difficult. However, it is possible to bring the challenges 
under some control by focusing on simple search strategies that cover multiple regions. For 
the displaced decays, just a few strategies are potentially sufficient. 


1. Search(es) for single vertex production, h —)• Gq + + ..., perhaps separated into 


37 



(a) h —» G 0 + + $t where the $ T is due to twin hadrons that decay invisibly and/or 
outside the detector. 

(b) h —> G 0+ + jet(s) where a promptly decaying twin hadron produces the jet(s). 

This type of search may only be feasible when requiring the presence of associated 
objects that may accompany the h, such as a lepton or a pair of vector boson fusion 
(VBF) jets. 

2. Search for exclusive di-vertex production production: as in h —> G 0 +G 0 +. 

3. Search for inclusive di-vertex production: h — >> 3 twin hadrons, of which at least two 
decay visibly and displaced. 

In contrast to exclusive di-vertex production, here the pair of observed twin hadrons 
generally have invariant mass below m,h, and need not be back-to-back in the h rest 
frame. 

To a limited degree, each of the three search strategies has been explored by ATLAS m 
5UU52H55] . CMS PHESHHE] and/or LHCb [3T]- However, due often to trigger limitations 
or analysis gaps, even the most sensitive of these searches do not yet put significant bounds 
on this model, and a broader and deeper program of searches is needed in Run II. We will 
discuss these issues further in Appendix [Cj 

We note also that although our minimal model has specific relationships between masses, 
lifetimes and production mechanisms, these relationships will not necessarily hold in other 
Twin Higgs and Twin Higgs-like models. It is therefore preferable that the above searches 
be carried out with the masses and lifetimes of the long-lived particles, and characteristics 
of the $ T (if any), treated as free parameters. In Appendix [c] we suggest benchmark models 
for these searches and consider some important triggering and analysis issues. 

B. New Higgs Decays Without Displaced Vertices 

When the G 0 + is heavy, its decays are prompt. Prompt non-SM decays of the Higgs 
such as h — y G 0 +G 0 + —^ (&&)(&&), ), or h — y Go+ T t (^) -f~ (j'+t ) T 

$ T , (/i + /i _ ) + fj T -, etc., were examined in a recent overview of non-SM Higgs decays [59]. 
Not all cases have yet been investigated for a 13-14 TeV machine, but it appears that a 
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non-SM branching fraction of order ICO 3 is far too small for discovery of such modes at 
LHC, though perhaps an e + e _ machine producing at least 10 5 Higgs bosons could find it. 
However, enhancement of the twin hadron branching fraction to ~ 10%, through h —> bb 


or non-perturbative effects, as discussed in Section |VE[ could bring these processes within 
reach of the LHC. 

At the other extreme, one can discuss completely invisible decays. These can dominate 
in region C, and can become important in regions A and B if the Go+ has an extremely long 
lifetime, or in region D if the G 0 + mixes substantially with invisible \ states. Other twists 
on the model (such as the presence of a massless twin hypercharge boson with small kinetic 
mixing, see Appendix [E]) can cause all twin hadrons to decay to invisible hidden particles. 
Detecting an invisible decay rate much smaller than 10% is very difficult at the LHC, and 
thus can only be done if there is significant enhancement by the h —> bb process. Again, an 
e + e _ collider would do much better. 


C. Heavy Higgs Decays 

Not only the h but also the heavy twin Higgs h may serve as a portal, if the twin SU (2) 
is linearly realized. The existence of a second perturbative Higgs, while not guaranteed, is 
favored by precision electroweak data [60]. If present it provides additional opportunities to 
uncover the twin mechanism. The mass of this second Higgs is ~ a/2A /, and so typically 
lies around the TeV scale. It possesses a v/ f -suppressed coupling to top quarks and thus 
is produced through gluon fusion, albeit with a (u//) 2 -suppressed cross-section compared 
to an SM Higgs of equivalent mass — falling at •> f~s = 14 TeV from 1000 fb at a mass of 
300 GeV to 10 fb for a mass of 1 TeV. 

On the lower end of this mass range, decays into WW,ZZ, and hh dominate, with 
branching ratios roughly proportional to 2 : 1 : 1. Once h —> WW, ZZ is kinematically 
allowed, these processes become at most comparable to WW, ZZ, and hh; although couplings 
to the SM bosons are suppressed by v/f due to mixing, the longitudinal coupling scales as 
m?/u 2 and entirely compensates. This is inevitable, since in the Twin Higgs mechanism the 
h and longitudinal W, Z, W, Z are all Goldstone modes under the same symmetry breaking. 

The decays h —s- hh,ZZ are common in many BSM models, and searches for these 
promising signals are already underway [6Ti 104] . The h will appear as a resonance with 
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width suppressed by (v/f) 2 compared to a SM Higgs of the same mass, times | (|) to 
account for the channel h —> hh (and channels h —y WW, ZZ, if kinematically allowed.) 
Observation of this resonance at equal rates in h —>• ZZ and h —y hh, and measurement of 
its width, could therefore allow for a test of the model. In Run II, CMS expects to exclude 
cr x Br (pp —* h —* ZZ) > 10 (4) fb for a heavy Higgs of 500 GeV (1 TeV) with 3000 ftr 1 
at a fs = 14 TeV, or to discover a x Br (pp —> h —> ZZ) > 30(10) fb at 5 a [65 j - This 
corresponds to an exclusion reach of 900 GeV or a discovery reach of 750 GeV for a h with 
/ = 3u. 

The h may also give rise to twin hadrons, either via h —> tt (followed by t —>■ bru) or 
h —* ZZ (with Z —y bb). Although the rate is less than the perturbative irreducible rate 
for h —y twin hadrons, it can easily happen that the trigger efficiency for h decays is low 
while that for h decays is high, so that h decays may provide an easier signal. Moreover 
the rate to produce twin hadrons via Z —y bb is larger by a factor of 10 or more than 
h —> ZZ —> a + a-l + d- and h —> hh —> 6677 , the cleanest h —> SM processes. Decays of 

h may often produce displaced glueballs, including ones too heavy to appear in h decays 
and decaying via an off-shell h [2D I23j. Even if all twin hadron decays are prompt, the 
events might be observable at LHC if the rate, multiplicity, and total (or missing) energy 
are sufficient [201 [ 66 ]. 

On the other hand, if the twin SU ( 2 ) is non-linearly realized and the h is as indistinct as 
the cr of QCD, then twin hadrons can be produced, albeit with small LHC rates, through 
enhancements of gg —» ZZ and VBF production of Z pairs. 

D. Precision Higgs Measurements 

Here we consider the role of the canonical signature of the Twin Higgs [Tj, namely 
0(v 2 /f 2 ) changes in Higgs couplings due to the misalignment between the electroweak vac¬ 
uum expectation value and the pseudo-Goldstone Higgs. This leads to a suppression of all 

2 

Higgs couplings by an amount 1 — fp relative to Standard Model predictions. There may 

also be a shift in branching ratios due to the additional partial width of decays into the 

twin sector, but as we have seen this can be much less than 10%. Assuming that Br (h —y 

twin sector) 10%, then the sole effect of the twin sector on SM Higgs measurements is 

2 

a reduction in all production rates, relative to SM predictions, by a factor of 1 — To 
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FIG. 8: Current bounds on v/f and the Higgs branching ratio into the twin sector from a combined 
fit to Higgs coupling measurements. Solid, dashed, and dotted black lines denote the 1-, 2-, and 3 -a 
bounds (defined as Ay; 2 = 2.30,6.18,11.83) due to ATLAS and CMS Higgs coupling measurements. 
The grey lines correspond to the perturbative calculation of the Higgs branching ratio into the twin 
sector as a function of v/f for yb/yb = 0,1,2; as discussed in the text, the actual branching ratio 
may differ significantly from the perturbative result for a given value of i)b- 

evaluate the impact of current Higgs coupling measurements on v/f, we have performed a 
combined fit of the most recent ATLAS and CMS Higgs measurements [(171 f7ol using the 
profile likelihood method m- The resulting bounds on v/f are shown in Fig. [8]as a function 
of v/f and the Higgs branching ratio into the twin sector. 9 We also show contours corre¬ 
sponding to the perturbative calculation of Br (h —y twin sector) as a function of v/f for 
Vb/Vb = 0,1,2. As discussed in Appendix [Bj the complications of bottomonium production 
suggest that the actual branching ratio is potentially much smaller than the perturbative 
value for sufficiently large Vb/Ubi while the irreducible rate for glueball production applies 

9 This fit does not include implicit precision electroweak bounds from infrared contributions to S and T. 
However, as we will discuss more in Appendix|Dj in contrast to composite Higgs models where the infrared 
contribution is cut off by m p ~ few TeV and provides the strongest constraint on coupling deviations [60], 
here the infrared contribution is cut off by the mass of the heavy Higgs. For m^ < TeV these corrections 
to S and T are comfortably compatible with current precision electroweak bounds and do not strongly 
influence the coupling fit. 
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for yb/yb = 0. Current measurements of Higgs couplings place a bound on v/f consistent 
with the benchmark value f — 3v considered here. 

For / ~ 3v the shift may be detected definitively before the end of the LHC, but not 
soon — certainly not within Run II. Future projections are a somewhat delicate matter, 
as measurements in certain channels will become systematics-limited and naive combina¬ 
tions neglecting correlated systematics are no longer appropriate. However, the collabora¬ 
tions quote appropriate coupling projections taking these effects into account. For example, 
ATLAS projects sensitivity to a uniform scaling factor k of 3.2% (2.5%) assuming cur¬ 
rent (no) theory systematics with 300 fb _1 M, which corresponds to a 95% CL bound of 
v/f < 0.35 (0.31) assuming the invisible branching ratio gives a sub-leading correction to 
the observed Higgs couplings. The remaining high-luminosity run should improve sensitiv¬ 
ity to 2.5% (1.6%) assuming current (no) theory systematics with 3000 fb -1 [77], leading to 
bounds of v/f < 0.31 (0.25). It is therefore unlikely that Higgs coupling measurements at 
the LHC could be used to substantially constrain the parameter space of Twin Higgs models, 
though if the twin mechanism is operative it may lead to 1 — 2cr deviations by the end of 
the high-luminosity run. (Direct limits on invisible decays are likely to fare even worse, with 
realistic projections suggesting 95% CL bounds will approach 8% with 3000 fb -1 [77], well 
short of the typical rate for decays into the hidden sector.) 

In any case, an overall coupling reduction is purely a sign of mixing and is non-diagnostic; 
while observation of this reduction would be revolutionary, no unique interpretation could 
be assigned to it. Only discovery of the twin hadrons and/or discovery and study of the 
mostly-twin heavy Higgs would allow for a clear interpretation of such a measurement. 

VII. CONCLUSIONS 

The “Twin Higgs” mechanism provides an existence proof for the unsettling possibil¬ 
ity that the solution to the hierarchy problem involves a sector of particles that carry no 
Standard Model quantum numbers, and are therefore difficult to produce at the LHC. The 
existence of dark matter already motivates us to consider hidden sectors, and it is important 
that the possibility of hidden naturalness be thoroughly considered. 

Fortunately, hidden sectors are often not as hidden as they first appear. As we have seen, 
a Fraternal Twin Higgs model, one whose hidden sector is not a precise twin of the Standard 
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Model, but which contains the minimal ingredients for the theory to address the hierarchy, 
naturally leads to Hidden Valley phenomenology. Specifically, the Higgs sector requires a 
twin top, which in turn favors twin color; and the lack of light twin quarks in this minimal 
model then leads to confinement and twin hadrons that include twin glueballs, the lightest 
of which necessarily mixes with the Standard Model-like Higgs, and twin quarkonium, whose 
tower of 0 ++ states has analogous mixing with the Higgs. We have shown that the resulting 
mixing can often cause these glueballs, and in some cases the lightest quarkonium states, 
to decay on an observable timescale, leading to a new source of visible non-Standard-Model 
Higgs decays. The branching fraction for these decays is at worst 10 -4 , typically at least 
1CF 3 , and possibly as large as allowed by Higgs measurements. For heavier twin hadrons 
with prompt decays, the visible final states — four heavy-flavor fermions, or two heavy-flavor 
fermions plus missing energy — are among those summarized in a recent overview of non-SM 
Higgs decays [59]. Searches for these final states and for invisible decays at LHC are possible 
if the branching fractions approach 10%; if much smaller, a lepton collider may be needed. 
But for moderately light glueballs (and perhaps quarkonia), these decays can produce a 
potentially spectacular signal of one or more highly displaced vertices, often accompanied 
by moderate amounts of missing energy from other twin hadrons that escape the detector. 
While several experimental searches for similar vertices have already taken place, a wider 
array of more powerful search strategies is required if the parameter space for this model is 
to be fully covered. 

The minimal Twin Higgs model we have presented is only mildly tuned, and no more 
unnatural than bottom-up modeling of Higgs compositeness or (natural) supersymmetry. 
This model also has no obvious cosmological problems, flavor problems, or other glaring 
issues. Thus we see no reason from the bottom up, given present knowledge, to treat 
Twin Higgs as less motivated than, say, composite Higgs models. Meanwhile, in addition 
to its cousin Folded Supersymmetry, the Twin Higgs model has recently been generalized, 
by recognizing it as an orbifold model [3B], 39]. Variants of these generalized models, by 
their very construction, share features with our minimal model, though they are different in 
details. This shows that the model space of sibling Higgses has not yet been fully explored, 
and should provide additional motivation for considering seriously and more generally the 
possibility of a hidden-sector solution to naturalness. 

On general grounds, we expect anything vaguely resembling a Twin Higgs model with a 
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hidden sector to require, as part of its solution to the naturalness problem, Higgs mixing via 
a “portal”-type interaction. This feature easily leads to additional Higgs-like resonances, 
new sources of missing energy, and exotic phenomenology of hidden-valley type, including 
non-SM Higgs decays to multi-body final states and/or displaced vertices. The challenges 
are that no individual model is required to produce any or all of these signals, and that 
production rates for these phenomena are not determined by known interactions (in contrast 
to gluino or stop production) and can be small. Among the most motivated places to search 
for new signals are in decays of known particles, whose production rates are known and large. 
New decays of the Higgs may not even be rare. In non-minima! models there may also be 
opportunities in rare decays of the Z. Searches for new phenomena generated by rarely 
produced heavy particles (such as the heavy Higgs in our model) must also be considered. 

Our field has tended to assume that the solution to the hierarchy problem lies in particles 
that resemble the ones that we know. While hidden sectors are often found in string theory 
vacua and required in models of supersymmetry breaking, their role has been limited to 
higher energy or purely gravitational interactions, leaving them, as far as the LHC is con¬ 
cerned, out of sight and out of mind. But with the possible scale of supersymmetry receding 
upwards, and with no sign of Higgs compositeness or of the colored top partners that were 
widely expected, the possibility of something more radical, such as a hidden sector around 
the weak scale that communicates with our sector through a portal, cannot be ignored. Our 
searches must move beyond the easier and more obvious lampposts, for the secrets of nature 
may lie hidden in the dark. 
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Appendix A: Quarkonium Mass Spectrum and Decays 


1. Spectrum of the Quarkonium States 


As explained in Section IV B , the most phenomenologically interesting quarkoninm states 
are the tower of 0 ++ states x, which can mix with the SM-like Higgs h and potentially 
decay back to the SM. We can make a crude but useful analysis of these states assuming 
they are governed by non-relativistic quantum mechanics; this approximation holds for low- 
lying x states in the regime rn- h A 3 . The effective quark-antiquark potential can be 
modeled as a combination of a (logarithmically corrected) Coulomb potential and a long¬ 
distance linear potential. However, x decay to SM particles is only of interest if x -»■ 
G 0 +G 0 + is kinematically forbidden, since otherwise the twin hadronic decay will dominate. 
This requires my < 2mo ~ I 3 . 6 A 3 . One may then check that (remembering a p-wave 
state is larger than an s-wave state) the Coulomb potential would imply an interquark 
distance of order A 3 1 for any relevant portion of the (mo, my) plane, showing the Coulomb 
approximation is very poor. This motivates using a purely linear potential, an approximation 
which improves for heavier states. 

Taking the potential as linear in r with a string tension a « 4Ay, the Schrodinger equation 
for the p-wave radial function is: 


1 / 1 d 2 d 2 

X - I n~ r — I 2 


2/i 


r* 


+ (err - E n ) 


R(r) = 0, 


(Al) 


r 2 dr dr 

with reduced mass ji — AA No exact solutions of this equation are known, but an excellent 
approximation may be obtained by neglecting the linear potential near r = 0 and the p-wave 
repulsive potential at large r, and then matching the two asymptotic oscillating solutions. 
Making this approximation we find the energy states 


/37r\ 2/3 A 4 / 3 
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n + - 
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2/3 


(A 2 ) 


and associated wave functions. This approximate solution is quite close to the exact numer¬ 
ical solution and may be used as a basis for more accurate estimates, though we have not 
done so here. In any case we believe that this formula correctly captures the parametric 
behavior of the energy levels and masses. 

In the Standard Model, highly excited bottomonium states above 10.56 GeV promptly 
decay to a pair of B mesons; a similar story applies for charmonium. But in this model there 
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are no light twin quarks, and correspondingly no twin B mesons, so the twin bottomonium 
states remain narrow until their masses are above 2m 0 or 2mg + m 0 . In most parts of 
parameter space these are very heavy scales, so that the towers of narrow states extend to 
much higher n than in the SM. 

For the physically interesting values of mo and mg, our approximation does not survive 
to high n, because relativistic effects become large at rather small n. Still, at small n 
our approximation should work for mg A 3 , including region A, most of the region B, and 
some portion of regions D and C. Full study of the spectrum with Coulombic and relativistic 
corrections is beyond our scope, but we expect that many of the results that we describe 
here will be parametrically valid (with order-one corrections) in the relativistic regime as 
well. 


2. Decays 


Here we explore decays of the light x states to the SM via mixing with the Higgs, and to 
a spin-one T state plus twin leptons l = f , v via a radiated off-shell Z. 

In models where the twin leptons are too heavy to allow \ —» T 11, the lowest-lying x 
state will dominantly decay to the SM, with a lifetime obtained from Eq. [43] 

'30 GeV\ 7 ( tyk \ 10/3 /750 GeV\ 4 


ct ~ 3 mm x ( - — W --- (A3) 

my J \2m 0 J V / ) 

for my > 10 GeV. Just as for glueball states (see Fig. [ 6 |, the lifetime jumps upward once 
my falls below the SM bb threshold. This formula will also apply for the n th x state (with 
an additional n~ 2 ^ suppression) if its twin strong and twin weak decays are kinematically 
forbidden. In this case, all the lowest-lying x states decay visibly even in region C, and 
displaced decays occur in low-mg portions of both regions C and D. 

Now suppose the twin leptons are light enough that y —>• TiJ? is kinematically allowed. 
(Note the T states will then always promptly annihilate into twin lepton pairs, which are 


invisible in this model, via an off-shell Z .) We saw in (44) that the decay x —)• Tgoes 


as mo — m -; 


While we could attempt to compute this splitting in the non-relativistic 


approximation, which would not be unreasonable for the higher n states, we know from SM 
charmonium that this splitting is strongly affected by spin-related effects that push my up 
and my down. Calculation of this splitting is therefore beyond the scope of our work. 
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FIG. 9: Features of x decay; we emphasize that these results have very substantial uncertainties 
due to the limitations of our estimation methods. Left: Branching ratios of \ into the SM along 
the BD-boundary line m.G 0+ = mg, assuming twin neutrinos are massless. Right: Decay length of 
X for mc 0+ = mg, where the lower red (upper blue) curve refers to massless (heavy) twin neutrinos. 


However, we may use a trick to get a good estimate in a phenomenologically important 
part of parameter space, namely the part of region D (where both bottomoninm and glneballs 
may be produced, often in the same event) that lies close to region B. This “BD-boundary” 
lies at niG 0+ ~ mg ~ |my, a mass ratio very similar to the SM charmonium system, for 
which lattice studies give a SM glueball mass of 1.5 — 1.7 GeV [03, 78, [79], and a charm 
quark mass m c ^ 1.5 GeV. Therefore we can take results from charmonium data [80]. 


m 


XcO 


= 3.415 GeV, mj /</, = 3.097 GeV, 


(A4) 


and rescale all mass scales by a common factor to obtain a realistic value for my — my near 
the BD-boundary. 

Using this scaling argument and the formulas of Section VD we see from Fig. [9] that 
displaced decays are possible for my somewhat below 30 GeV, with lifetimes potentially 
reaching 10 cm; otherwise the decays are prompt. The branching ratio of the lightest x f° 
SM states is larger than ~ 10% (larger than ~ 1%) for 2mg ~ 25 GeV (10 GeV). We also 
remind the reader that in this part of region D one expects the Higgs branching rate to 
twin bottomoninm, from the h —)■ bb process, to be well enhanced above the irreducible rate. 
Consequently even a rather small branching fraction can be interesting. 

However, we must emphasize that our estimates are highly uncertain due to our crude 


methods, which include the fact that we have not computed the numerical coefficient in (44). 
One should therefore only conclude that there is good reason to believe that the lowest x 
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may be observable near the BD boundary, and possibly displaced, even if twin leptons are 
light. Excited x states will decay invisibly. 

As we increase m 0 , moving away from the BD-boundary across region D and into C, we 
cannot reliably compute the m 0 /mg dependence of the x ~ T splitting and corresponding 
twin weak decay rates. However, both our approximations and reasoning from QCD data 
suggest the mass splitting rapidly increases with m o, at least as fast as A 3 oc mo itself. We 
therefore expect that even the lowest x state decays invisibly with a branching fraction near 
unity, across the right side of region D and all of region C. 

Finally let us note that we are not aware of any other bottomonium states that have 
observable decays. 

Appendix B: Twin Hadron Production in More Detail 

Twin hadrons are produced when h —y gg or bb , but twin hadronization is complicated and 
poorly understood. Despite many years of experimental study of QCD hadron production, 
theoretical understanding of hadronization is still limited. Moreover a twin sector with no 
light twin quarks has a very different spectrum and dynamics from SM QCD. 

Perturbation theory applies for inclusive calculations, such as Br {h — y twin hadrons), 
when colored particle production occurs at distances short compared to Ag 1 , so that a s is 
small. (A similar example is the hadronic branching fraction of the Z.) However, for large 
mo or m^, twin hadron production is analogous to production of QCD hadrons within or just 
above the QCD resonance region, where non-perturbative effects can be important. Decays 
of the h can be enhanced relative to perturbative estimates if lies close to a narrow 
excited 0 ++ glueball or quarkonium resonance, and suppressed if it lies between two such 
resonances. 


1. Nonperturbative effects in twin gluon production 


First consider h —y gg for mg > m^/2. For small m 0 , the glueball resonances near 
125 GeV are likely to be both numerous and wide, blending into a continuum. In this case 


the irreducible rate is given by the perturbative result (45) for h —>■ gg 


For larger mo, resonance effects are potentially important. Based on large- N c counting 
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and experience from the p meson — Y p ~ m p j 5 is of order 1/N C = 1/3 — we expect that 
excited 0 ++ glueball resonances have Y/m ~ 1/Af? = 1/9, perhaps with an additional minor 
suppression factor. There is also likely a phase-space effect increasing the width for higher 
excitations, but we ignore this for our conservative estimate. Meanwhile the spacing between 
0++ glueball resonances is likely to be smaller than mo (as is the case for p mesons, among 
others) but it is surely no larger; lattice evidence for the second 0 ++ glueball [42] supports 
this. The maximum suppression between two resonances of width T and spacing A, relative 
to a perturbative calculation, is of order T/A > (m 0 /N^)/m 0 = 1/9. Thus we expect the 
worst possible non-perturbative suppression factor is about 0.1, bringing the worst-case rate 
for glueball production down no further than 10“ 4 . 


2. Nonperturbative effects in twin bottom production 


If the Higgs can decay to b quarks and thus to [ bb] states, the rate for twin hadron produc¬ 
tion is often enhanced. As before, the [bb] production rate is given by the perturbative rate 
for h —)■ bb for sufficiently small mg and mo- As discussed in Section 


YE 


the perturbative 

branching fraction for h —» [bb] + X grows as y/, contradicting existing Higgs measurements 
for y b > 1.25 y b . 

But at high y b , the perturbative rate often gives the wrong answer. The widths of the 
excited x states with mass ~ rn b may be very small compared to their mass splittings, due 
to kinematic constraints. For instance, the x widths are tiny if rrih < mfj + mo, 2m,), and 
may be quite narrow until rn b — m p m o. If m/, lies between two resonances, then there 
can be a strong non-perturbative suppression compared to the perturbative rate. 

We may make an estimate of the maximal suppression factor as we did for glueballs for 
large my and small m 0 . Here the annihilation decay x ~> 99 is perturbative and depends on 
-R'(O), which we have estimated as described in Appendix[Aj this gives us T if radiative decays 
of the x are sufficiently suppressed by small twin bottomonium mass splittings. Similarly 


we may obtain A for the higher x states from (A2). We find a suppression factor of order 


«3(^K) 3/2 , 


(Bl) 


where a% arises from the annihilation rate and should be evaluated at the scale mg, and 
where E n is the excitation energy of the relevant x state. Since the relevant x state has a 
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mass ~ rrihi E n ~ rrih — 2 m b , as long as n is not too small. This suppression factor can be 
of order 0.01-0.1. A suppression factor of (i/ft/1.25) 2 ~ (2m^/30 GeV) 2 would be enough to 
make a perturbatively excluded region allowed, so we conclude that some regions with m b 
somewhat below 2, specifically those where mj t lies between two excited bottomonium 
states, are probably not excluded by data. 

For larger mo there can be even more suppression because the rate x ~* 99 can be mod¬ 
ulated by excited glueballs for the same reason as h —> gg is modulated by excited glueball 
resonances. We therefore expect that the non-perturbatively allowed portion of region A 
extends to lower m b than that of region B. However, we have no reliable computational 
methods in this regime. 

Appendix C: Search Strategies 

Here we briefly discuss some triggering and analysis issues with regard to the search for 
the displaced decays of long-lived twin hadrons. For brevity, we refer to the visible decay 
products of a long-lived particle as a “displaced vertex” (DV) even when it occurs in regions 
where tracks are not actually reconstructed. We also refer collectively to the decay products 
of a W or Z in Wh and Zh production, and to the vector boson fusion (VBF) jets in VBF 
Higgs production, as “associated objects” (AO). 

1. Comments on Triggering 

Triggering on Higgs decays to long-lived neutral particles requires three classes of trigger 
strategies, which respectively focus on 

• the presence of one or two DVs; this method is largely independent of how h is produced 
and is sensitive to non-Higgs production of the DVs. Displaced decay objects can 
include jets with displaced tracks, trackless jets (possibly including a muon), narrow 
trackless jets with little electromagnetic calorimeter (ECAL) deposition, and unusual 
clusters of hits or tracks in the muon system. 

• the presence of AOs that accompany the Higgs, including VBF jets or daughters of a 
W or Z (leptons, neutrinos, jets); this method is relatively independent of the details 
of the h decay and can therefore be used for any exotic h decay mode. 
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• the presence of both; for instance, in an event with VBF jets along with a trackless 
jet. Requiring both may be used to lower p T thresholds on the trigger objects, or to 
access DVs that would be unusable on their own. However, this powerful method is 
specifically optimized for Higgs decays to long-lived particles. 


Triggers of the first type were used at ATLAS |26l - f28l 1501 f53Tf55l [ST ] and CMS 
in Run I; lepton and $ T triggers are standard, while a VBF trigger was used in 2012 parked 
data at CMS [83j; and no triggers of the third type were used in Run I, to our knowledge. 

Depending on the lifetimes of the long-lived particles and on their masses and other 
kinematics, any one of the three approaches to triggering can work for any of the three search 
strategies outlined in Section |Y1 A| We leave the appropriate studies to our experimental 
colleagues. 


Benchmark Models 


Next we point out possible benchmarks that may be used as straw-man targets. Although 
the parameters of the benchmarks are correlated in the Fraternal Twin Higgs — for instance 


the mass and lifetime of the Gq+ are highly correlated, Eq. (40) — it is important to account 


for the many uncertainties in twin hadronization, and even more important, to retain model- 
independence. Therefore, notwithstanding the particular features of the Fraternal Twin 
Higgs, it seems best to ignore these correlations and study the benchmark models across 
their entire parameter space. 

In Section [VTA] we suggested three possible search strategies, and we now discuss suitable 


benchmarks appropriate for each of them. To keep things simple, we discuss them in the 
context of phenomena that occur in region A. However, with a little thought the reader may 
verify that the same benchmarks would be useful for these and other processes that occur 
in other regions of parameter space. 


a. Single Displaced Vertex Search 

To look for a single hadronic DV from h —y Gq + + ... (and similar decays) is challenging, 
because of large, difficult-to-measure backgrounds. A CMS search [59] that required only 
one DV was not able to set limits on decays of a 125 GeV particle. Relevant ATLAS searches 
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used single DV events to obtain background estimates on double DV events. However, such a 
search is necessary. When ctq 10 m (as for m 0 10 GeV in region A), or if hadronization 
assures that particles making DVs are rare among twin hadrons, the number of DVs rarely 
exceeds one per event. 

To date no search for a DV has exploited the AOs, the VBF jets and/or lepton(s), that 
sometimes accompany the Higgs. To obtain background estimates, it may be enough to 
measure the rates for events with neither an AO nor a DV, and with either an AO or a DV; 
then if the AO and DV are uncorrelated one may predict the rate to have both an AO and 
a DV. 

Also, searches for a single DV can demand a second object from the twin hadrons that 
do not produce a DV. This object could be $ T (relevant for very long-lived G*o+) or prompt 
jets, possibly 6 -tagged (relevant for h —* Gq + G' +0i where G ' +0 —* bb promptly.) 

In this context we suggest the following benchmark models (in which all X t of mass m l 
decay to SM final states with the same branching fractions as a Higgs of mass m,, unless 
otherwise noted). These benchmarks do not cover all the kinematic possibilities but will 
serve as a useful initial target. 

1. h — > X 1 X 2 , where X\ has a long lifetime r±, and X 2 —> $?■ 

2. h — > X 1 X 2 , where Ad has a long lifetime T\ and now X 2 — > SM with a prompt decay. 

It is important to consider a range of masses for m\ and m 2 . For the first benchmark, m\ 7 ^ 
m 2 helps account for the possibility of h — > many glueballs, of which only one gives a DV 
and the rest are unseen; then X 2 represents a system of > 2 glueballs, with a large invariant 
mass. The signal of the second benchmark may arise from decays such as h —)■ Gq + G' 0+ or 
h —>■ Go+X, which generally have m\ 7 ^ jn 2 . 

b. Exclusive Double Displaced Vertex Search 

In the heart of region A (m 0 ~ 20 — 40 GeV), h —> G 0+ G 0+ often gives two DV, roughly 
back-to-back and with invariant mass ~ Several searches for this final state have been 
undertaken at ATLAS mmm and at LHCb ph Also interesting is h —>■ Go+G' +0 ; here 
one has particles of two different lifetimes and masses. 

In this context we suggest benchmark models of the form 
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1. h —> XX, where A" has mass nix and lifetime tx- 


2 . h —>■ XiX 2 , where now the A* have different masses mi < m 2 and lifetimes r ? > 

In onr minimal model it is sufficient to consider m i < m 2 < 2 m i since otherwise A " 2 —>■ Ad Ad 
decays occur. This will not be true in all models however. 

c. Inclusive Double Displaced Vertex Search 

For mo ~ 10—25 GeV, with lifetimes still within the detector, the probability of producing 
> 2 glueballs per event is larger, making it more common to have two Go+ DVs along with 
one or more invisible glueballs. One key difference from the exclusive double DV search is 
that the two DVs may not be well-separated in the lab frame, and may even tend to be 
found in the same angular region of the detector. Also, the invariant mass of the two (or 
more) DVs may be well below m/,.. 

A suitable benchmark highlighting this difference would be 

• h —> A 1 A 1 A 2 , where the masses are chosen freely (but typically m 2 > mi), Ad 
decays as usual with lifetime T \, and X 2 —>■ $ T . Here the two DVs will have a broad 
distribution in angular separation and in the momentum that they carry in the h rest 
frame. 

The one serious caveat is the possibility (perhaps remote in the our minimal model, 
but not necessarily in other models) that the DVs are often clustered pA], 16| [20, 52) [83J. 
A special benchmark model may be needed if this clustering is sufficiently common that 
isolation requirements on DV candidates typically fail. We suggest 

• h —> X 0 X 2 , followed promptly by A 0 — > Ad Ad, Ad decays as usual with lifetime T\, 
and X 2 —>- $ T . By choosing 2mi < mo < m 2 <C m/,, one assures that A " 0 is relativistic 
and the two Ad decays are correlated in angle. 

Appendix D: Precision Electroweak 

Here we briefly consider precision electroweak constraints on the Twin Higgs mechanism. 
In any theory where coupling deviations of the SAI-like Higgs are due to mixing with heavier 
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FIG. 10: Close-up of the one- and two-sigma S, T ellipses marginalized over U with contributions 
from both the SM-like Higgs and the heavy twin Higgs h. The small (large) blue “x” shows the 
S, T prediction including a heavy twin Higgs of mass m; = 750 GeV (1.5 TeV) for / = 3v, while 
the small (large) red “+” shows the S, T prediction including a heavy twin Higgs of mass mt = 750 
GeV (1.5 TeV) for / = 5v. The black dot denotes the contribution of a purely Standard Model 
Higgs of mass rrih = 125 GeV. 

states there are infrared contributions to the S and T parameters whose coefficients depend 
on the reduced coupling of the SM-like Higgs to electroweak gauge bosons. The logarithmic 
contributions to S and T scale like log(m^/mz) + (l—« 2 ) log(A/m/J, where a = ghvv/9h SM vv 
and A is the cutoff scale associated with heavy states. In composite Higgs models the 
coupling deviations come from mixing with heavy resonances, which must lie in the multi- 
TeV range due to tree-level contributions to S and T. The cutoff on infrared contributions 
is then of order A = A'kv j\J\ — a 2 ~ 47r/, leading to extremely strong constraints on a [GO] . 

In the Twin Higgs mechanism, the 0(y/f) deviations in Higgs couplings come from 
mixing with the heavy twin Higgs and the cutoff is precisely A = rri- h . If the additional twin 
Higgs is nonperturbatively heavy (i.e. rn- h ~ 47r/), the situation is identical to composite 
Higgs models, but in general the heavy twin Higgs can be much lighter thanks to a custodial 
symmetry; there are no tree-level contributions to S and T associated with h. The complete 
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contributions to S and T from the SM-like Higgs and heavy Higgs take the form 
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relative to the precision electroweak contributions of a theory with a SM Higgs of mass mh- 
The current GFitter bound gives S = 0.05 ± 0.11 and T = 0.09 ± 0.13 for a reference Higgs 
mass of mh = 125 GeV |85j . We translate this into a fit by marginalizing over U to create 
the S, T ellipse for free U and plot the spread of IR S, T contributions as a function of m- h 


and / in Fig. 10 Remaining within the la (2cr) error for f = 3v requires rn- h < 600 GeV 
(5 TeV). In contrast to composite Higgs models, precision electroweak constraints do not 
substantially constrain Higgs coupling deviations in the Twin Higgs as long as the heavy 
Higgs is in the TeV range. 


Appendix E: The Effect of Twin Hypercharge 

In our minimal model, for reasons described in Section [TTJ we have assumed twin hyper¬ 
charge is absent. What if it is present? Then the most natural values of the mass for the 
twin photon 7 are (1) ~ A, in which case nothing changes from our present discussion, or 
( 2 ) zero. 

A key issue for twin photons is their kinetic mixing with the SM photon [HSj- 

In our model the mixing parameter n is not generated at low-loop order at scales below 
A, so its size is determined by its value at A, which is itself determined by physics in the 
ultraviolet that we do not specify. It is thus a free parameter of the model. 

If rriA, = 0, we define the twin photon as the massless particle that does not couple to 
the visible sector, and the f and b acquire electric charges e = ne/e and e/3, respectively, 
in units of e. The twin photon is stable and invisible. All twin hadrons may decay via 
radiation or annihilation to on-shell or off-shell twin photons — twin bottomonium through 
tree-level graphs and twin glueballs through a loop of twin bottom quarks. The rates for 
these decays far exceed rates for decays through the off-shell Higgs, and so most decays to 
the hidden sector are invisible. 

However, any twin photon may be replaced by a SM photon at the cost of e 2 e 2 in branch¬ 
ing fraction. Constraints on £ for a massless (or nearly massless) twin photon arise from 
e + e _ production of the fractionally charged twin f and b fermions. Since most twin sector 
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production events are invisible, one can obtain a model-independent limit on e that depends 
only on whether twin f and b are kinematically accessible, for instance in Z decays. Such 
limits are weak, however, of order e < 1CT 1 [87]. But typically the true limits are stronger, 
because the fraction ~ £ 2 of events that produce a visible signal could have been seen at 
LEP. The detailed limits are far too complex for us to work out here. 

If e is not small, lying in the 10 ~ 3 — 10 -1 range, h — > 7 + $ T and rare h —» 77 + $ T 
decays may occur, especially if the twin hadron production rate is large (regions B, C, and D 
with m~ b not too small). Moreover, annihilation of spin-one states via an off-shell SM photon 
may produce a visible final state, e.g. T —» 7 * —>• £ + £~. This prompt resonant decay will 
be subleading by ~ e 2 compared to the invisible decay T —>■ 7 * —$■ tt as long as m? < m b ] 
otherwise, it competes favorably with highly-suppressed twin sector decays, and may even 
dominate. This type of decay is probably common in regions C and D. In region B, where 
bottomonium decays to glueballs, this decay can perhaps still occur for the 1 glueball, 
but this glueball is heavy and probably quite rare. 

In summary, for rriy = 0 and at small e, all twin-sector decays are invisible, but partly- 
visible decays arise at large e. In region A, where production is small, + decays may be 
the best bet but are very rare. In region B, the larger production rate offers the possibility 
of a large 7 + $ T rate and a small but important 77 + $ T rate. And in regions C and D, 
the final state + p T , where the /i + /i~ form the lowest T resonance, may be either 

rare or common. (Higher T resonances will often decay radiatively.) Thus searches for both 


invisible and partly visible decays of the Higgs are highly motivated; see Section VIB above. 

Meanwhile, a twin photon with nn.y -C A typically introduces a new naturalness problem 
if Higgsed, since one more Higgs field S is required and its expectation value and mass 
must be -C A. One can get around this if the hypercharge coupling g± is very small, with 
Pi A < rrih/2. Alternately, the twin photon may be rendered massive by the Stueckelberg 
mechanism, which is technically natural but entails a parametrically new scale if my -C A. 
In either case, we view this possibility as lying beyond the minimal Twin Higgs story. 

For 2 m e < m^ < mh , the twin photon can play a major role in LHC signals. In contrast to 
the massless case, here we define the massless SM photon as decoupled from twin particles; 
the twin photon is identified as the mass eigenstate, which couples to SM fields with strength 
~ ee. Assuming the twin photon is light enough that all twin hadrons can decay to it, twin 
hadron production generically (and not just for 0 ++ hadrons) leads to high multiplicity 
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(> 4) twin photon final states. Through kinetic mixing with the SM photon, the twin 
photons then decay to lepton pairs and jet pairs (or simply hadron pairs, such as K + K~ 
and 7r + 7r _ , if ~ few GeV or below.) The decays of the twin photon may be prompt 
or displaced depending on the value of e, on which there are a number of constraints; for 
a recent summary, see [88]. If prompt, however, limits on multilepton final states (for 
rtu, 1 GeV) [891-192] and on prompt “lepton-jets” (for rru ( < 1 GeV) [931198] are strong 
constraints. Searches for significantly displaced lepton pairs with masses well above 1 GeV 
[55 J and < 1 GeV [52] [5H, 55] are sensitive to this signal but have by no means excluded it. 
Moreover, any such constraints become much weaker once ct.^ ~ 10 m. Note also 

that if some twin hadrons are too light to decay to twin photons, a mixture of these signals 
and those described in the main text is possible. 


[1] Z. Chacko, H.-S. Goh, and R. Harnik, The Twin Higgs: Natural electroweak breaking from 
mirror symmetry, Phys.Rev.Lett. 96 (2006) 231802, hep-ph/0506256|. 

[2] R. Barbieri, T. Gregoire, and L. J. Hall, Mirror world at the large hadron collider, 

hep-ph/0509242. 

[3] Z. Chacko, Y. Nomura, M. Papucci, and G. Perez, Natural little hierarchy from a partially 
goldstone twin Higgs , JHEP 0601 (2006) 126, |hep-ph/0510273 . 

[4] G. Burdman, Z. Chacko, H.-S. Goh, and R. Harnik, Folded supersymmetry and the LEP 
paradox, JHEP 0702 (2007) 009, |hep-ph/0609152 . 

[5] G. Burdman, Z. Chacko, R. Harnik, L. de Lima, and C. B. Verhaaren, Colorless Top 
Partners, a 125 GeV Higgs, and the Limits on Naturalness, arXiv : 1411.3310. 

[6] S. Dimopoulos and G. Giudice, Naturalness constraints in supersymmetric theories with 
nonuniversal soft terms, Phys.Lett. B357 (1995) 573-578, |hep-ph/9507282|. 

[7] A. G. Cohen, D. Kaplan, and A. Nelson, The More minimal supersymmetric standard model, 
Phys.Lett. B388 (1996) 588-598, |hep-ph/9607394 . 

[8] R. Sundrum, SUSY Splits, But Then Returns, JHEP 1101 (2011) 062, arXiv : 0909.5430 . 

[9] R. Barbieri, E. Bertuzzo, M. Farina, P. Lodone, and D. Pappadopulo, A Non Standard 
Supersymmetric Spectrum, JHEP 1008 (2010) 024, arXiv : 1004.2256 1 . 

[10] M. Papucci, J. T. Ruderman, and A. Weiler, Natural SUSY Endures, JHEP 1209 (2012) 


57 



035, [arXiv: 1110.6926]. 

[11] C. Brust, A. Katz, S. Lawrence, and R. Sundrum, SUSY, the Third Generation and the 
LHC , JHEP 1203 (2012) 103, [arXiv: 1110.6670 . 

[12] J. A. Evans, Y. Kats, D. Shih, and M. J. Strassler, Toward Full LHC Coverage of Natural 
Supersymmetry, JHEP 1407 (2014) 101, |arXiv: 1310.57581. 

[13] J. Berger, M. Perelstein, M. Saelim, and P. Tanedo, The Same-Sign Dilepton Signature of 
RPV/MFV SUSY, JHEP 1304 (2013) 077, arXiv: 1302.2146], 

[14] Z. Han, A. Katz, M. Son, and B. Tweedie, Boosting Searches for Natural SUSY with RPV 
via Gluino Cascades, Phys.Rev. D87 (2013) 075003, arXiv: 1211.4025|. 

[15] M. J. Strassler and K. M. Zurek, Echoes of a hidden valley at hadron colliders, Phys.Lett. 
B651 (2007) 374-379, [hep-ph/0604261 . 

[16] M. J. Strassler, Possible effects of a hidden valley on supersymmetric phenomenology, 

hep-ph/0607160. 

[17] M. J. Strassler and K. M. Zurek, Discovering the Higgs through highly-displaced vertices, 
Phys.Lett. B661 (2008) 263-267, |hep-ph/0605193 . 

[18] T. Han, Z. Si, K. M. Zurek, and M. J. Strassler, Phenomenology of hidden valleys at hadron 
colliders, JHEP 0807 (2008) 008, [arXiv: 0712.2041], 

[19] M. J. Strassler, Why Unparticle Models with Mass Gaps are Examples of Hidden Valleys, 
arXiv:0801.0629. 

[20] M. J. Strassler, On the Phenomenology of Hidden Valleys with Heavy Flavor, 

arXiv:0806.2385. 

[21] J. E. Juknevich, D. Melnikov, and M. J. Strassler, A Pure-Glue Hidden Valley I. States and 
Decays, JHEP 0907 (2009) 055, arXiv : 0903.0883 1. 

[22] J. Kang and M. A. Luty, Macroscopic Strings and ’ Quirks’ at Colliders, JHEP 0911 (2009) 
065, [arXiv:0805.4642], 

[23] J. E. Juknevich, Pure-glue hidden valleys through the Higgs portal, JHEP 1008 (2010) 121, 
[arXiv: 0911.5616]. 

[24] L. M. Carpenter, D. E. Kaplan, and E.-J. Rhee, Reduced fine-tuning in supersymmetry with 
R-parity violation, Phys.Rev.Lett. 99 (2007) 211801, |hep-ph/0607204|. 

[25] ATLAS Collaboration, G. Aad et. al, Triggering on Long-Lived Neutral Particles in the 
ATLAS Detector, Tech. Rep. ATL-PHYS-PUB-2009-082, ATL-COM-PHYS-2009-233, 2009. 



[26] ATLAS Collaboration, G. Aad et. al., Triggers for displaced decays of long-lived neutral 
particles in the ATLAS detector, JINST 8 (2013) P07015, arXiv: 1305.22841. 

[27] ATLAS Collaboration, G. Aad et. al., Search for a light Higgs boson decaying to long-lived 
weakly-interacting particles in proton-proton collisions at yfs = 7 TeV with the ATLAS 
detector, Phys.Rev.Lett. 108 (2012) 251801, jarXiv: 1203.13031. 

[28] ATLAS Collaboration, G. Aad et. al, Search for pair produced long-lived neutral particles 
decaying in the ATLAS hadronic calorimeter in pp collisions at yfs = 8 TeV, Tech. Rep. 
ATLAS-CONF-2014-041, 2014. 

[29] CMS Collaboration, Search for long-lived neutral particles decaying to quark-antiquark pairs 
in proton-proton collisions at yfs = 8 TeV, arXiv: 1411.6530. 

[30] M. R. Buckley, V. Halyo, and P. Lujan, Don’t Miss the Displaced Higgs at the LHC Again, 

arXiv:1405.2082. 

[31] Y. Cui and B. Shuve, Probing Baryogenesis with Displaced Vertices at the LHC, 

arXiv:1409.6729. 

[32] H. Cai, H.-C. Cheng, and J. Terning, A Quirky Little Higgs Model, JHEP 0905 (2009) 045, 
[arXiv: 0812.0843]. 

[33] A. Falkowski, S. Pokorski, and M. Schmaltz, Twin SUSY, Phys.Rev. D74 (2006) 035003, 
hep-ph/0604066. 

[34] S. Chang, L. J. Hall, and N. Weiner, A Supersymmetric twin Higgs, Phys.Rev. D75 (2007) 
035009, |hep-ph/0604076 . 

[35] N. Craig and K. Howe, Doubling down on naturalness with a supersymmetric twin Higgs, 
JHEP 1403 (2014) 140, [arXiv: 1312.1341], 

[36] P. Batra and Z. Chacko, A Composite Twin Higgs Model, Phys.Rev. D79 (2009) 095012, 

[arXiv: 0811.0394], 

[37] M. Geller and O. Telem, A Holographic Twin Higgs Model, arXiv: 1411.2974. 

[38] N. Craig, S. Knapen, and P. Longhi, Neutral Naturalness from the Orbifold Higgs, 

arXiv:1410.6808, 

[39] N. Craig, S. Knapen, and P. Longhi, The Orbifold Higgs, arXiv : 1411.7393. 

[40] N. Arkani-Hamed, A. G. Cohen, and H. Georgi, Electroweak symmetry breaking from 
dimensional deconstruction, Phys.Lett. B513 (2001) 232-240, |hep-ph/0105239|. 

[41] N. Arkani-Hamed, A. G. Cohen, T. Gregoire, and J. G. Wacker, Phenomenology of 


59 




electroweak symmetry breaking from theory space, JHEP 0208 (2002) 020, |hep-ph/0202089 . 

[42] N. Arkani-Hamed, A. Cohen, E. Katz, A. Nelson, T. Gregoire, et. al., The Minimal moose 
for a little Higgs , JHEP 0208 (2002) 021, |hep-ph/0206020 . 

[43] M. Schmaltz and D. Tucker-Smith, Little Higgs review, Ann.Rev.Nucl.Part.Sci. 55 (2005) 
229-270, [hep-ph/0502182]. 

[44] C. J. Morningstar and M. J. Peardon, The Glueball spectrum from an anisotropic lattice 
study, Phys.Rev. D60 (1999) 034509, jhep-lat/9901004|. 

[45] Y. Chen, A. Alexandru, S. Dong, T. Draper, I. Horvath, et. al., Glueball spectrum and 
matrix elements on anisotropic lattices, Phys.Rev. D73 (2006) 014516, hep-lat/0510074|. 

[46] M. Gockeler, R. Horsley, A. Irving, D. Pleiter, P. Rakow, et. al., A Determination of the 
Lambda parameter from full lattice QCD, Phys.Rev. D73 (2006) 014513, |hep-ph/0502212[. 

[47] B. Lucini, M. Teper, and U. Wenger, Glueballs and k-strings in SU(N) gauge theories: 
Calculations with improved operators, JHEP 0406 (2004) 012, |hep-lat/0404008|. 

[48] A. Djouadi, J. Kalinowski, and M. Spira, HDECAY: A Program for Higgs boson decays in 
the standard model and its supersymmetric extension, Comput.Phys.Commun. 108 (1998) 
56-74, hep-ph/9704448]. 

[49] K. Juge, J. Kuti, and C. Morningstar, Bag picture of the excited QCD vacuum with static Q 
anti-Q source, Nucl.Phys.Proc.Suppl. 63 (1998) 543-545, jhep-lat/9709132 . 

[50] ATLAS Collaboration, G. Aad et. al., Search for displaced vertices arising from decays of 
new heavy particles in 7 TeV pp collisions at ATLAS, Phys.Lett. B707 (2012) 478-496, 

jarXiv: 1109.2242], 

[51] LHCb Collaboration, Search for Higgs-like bosons decaying into long-lived exotic particles, 
Tech. Rep. LHCb-CONF-2012-014, CERN-LHCb-CONF-2012-014, 2012. 

[52] ATLAS Collaboration, G. Aad et. al., Search for long-lived neutral particles decaying into 
lepton jets in proton-proton collisions at y/s = 8 TeV with the ATLAS detector, JHEP 1411 
(2014) 088, |arXiv: 1409.0746], 

[53] ATLAS Collaboration, G. Aad et. al, Search for long-lived, heavy particles in final states 
with a muon and multi-track displaced vertex in proton-proton collisions at yfs = 7 TeV with 
the ATLAS detector, Phys.Lett. B719 (2013) 280-298, jarXiv: 1210.745l|. 

[54] ATLAS Collaboration, G. Aad et. al., Search for displaced muonic lepton jets from light 
Higgs boson decay in proton-proton collisions at y/s = 7 TeV with the ATLAS detector, 


60 


Phys.Lett. B721 (2013) 32-50, jarXiv: 1210.0435). 

[55] ATLAS Collaboration, G. Aad et. al, Search for Displaced Muon Jets from light Higgs 
boson decay in proton-proton collisions at sqrt(s) = 7 Tev with the ATLAS detector , Tech. 
Rep. ATLAS-CONF-2012-089, ATLAS-COM-CONF-2012-112, 2012. 

[56] CMS Collaboration, V. Khachatryan et. al, Search for ”Displaced Supersymmetry'’ in 
events with an electron and a muon with large impact parameters , arXiv: 1409.4789, 

[57] CMS Collaboration, CMS, Search for long-lived particles decaying to final states that 
include dileptons, Tech. Rep. CMS-PAS-EXO-12-037, 2014. 

[58] CMS Collaboration, S. Chatrchyan et. al., Search in leptonic channels for heavy resonances 
decaying to long-lived neutral particles, JHEP 1302 (2013) 085, arXiv: 1211.2472|. 

[59] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, et. al., Exotic Decays of the 125 GeV Higgs 
Boson, Phys.Rev. D90 (2014) 075004, jarXiv: 1312.4992]. 

[60] J. Espinosa, C. Grojean, M. Muhlleitner, and M. Trott, First Glimpses at Higgs’ face, JHEP 
1212 (2012) 045, (arXiv: 1207.1717]. 

[61] CMS Collaboration, S. Chatrchyan et. al., Search for the resonant production of two Higgs 
bosons in the final state with two photons and two bottom quarks, . 

[62] ATLAS Collaboration, G. Aad et. al., Search For Higgs Boson Pair Production in the 77 bb 
Final State using pp Collision Data at y/s = 8 TeV from the ATLAS Detector, 

arXiv:1406.5053, 

[63] CMS Collaboration, S. Chatrchyan et. al., Measurement of the properties of a Higgs boson 
in the four-lepton final state, Phys.Rev. D89 (2014) 092007, |arXiv: 1312.53531. 

[64] ATLAS Collaboration, G. Aad et. al., Measurements of Higgs boson production and 
couplings in diboson final states with the ATLAS detector at the LHC, Phys.Lett. B726 
(2013) 88-119, jarXiv: 1307.1427], 

[65] CMS Collaboration, 2HDM Neutral Higgs Future Analysis Studies, Tech. Rep. 
CMS-PAS-FTR-13-024, CERN, Geneva, 2013. 

[66] J. A. Evans, A Swarm of Bs, JHEP 1408 (2014) 073, arXiv: 1402.4481]. 

[67] ATLAS Collaboration, G. Aad et. al., Search for the bb decay of the Standard Model Higgs 
boson in associated ( W/Z)H production with the ATLAS detector, arXiv : 1409.6212, 

[68] ATLAS Collaboration, Search for the Standard Model Higgs boson produced in association 
with top quarks and decaying to bb in pp collisions at yfs = 8 TeV with the ATLAS detector 


61 




at the LHC, Tech. Rep. ATLAS-CONF-2014-011, CERN, Geneva, Mar, 2014. 

[69] ATLAS Collaboration, Evidence for Higgs boson Yukawa couplings in the H rr decay 
mode with the ATLAS detector, Tech. Rep. ATLAS-CONF-2014-061, CERN, Geneva, Oct, 
2014. 

[70] ATLAS Collaboration, Observation and measurement of Higgs boson decays to WW* with 
ATLAS at the LHC, Tech. Rep. ATLAS-CONF-2014-060, CERN, Geneva, Oct, 2014. 

[71] ATLAS Collaboration, G. Aad et. al, Measurements of Higgs boson production and 
couplings in the four-lepton channel in pp collisions at center-of-mass energies of 7 and 8 
TeV with the ATLAS detector, arXiv: 1408.5191. 

[72] ATLAS Collaboration, G. Aad et. al., Measurement of Higgs boson production in the 
diphoton decay channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the 
ATLAS detector, arXiv: 1408.7084. 

[73] CMS Collaboration, Precise determination of the mass of the Higgs boson and studies of the 
compatibility of its couplings with the standard model, Tech. Rep. CMS-PAS-HIG-14-009, 
CERN, Geneva, 2014. 

[74] CMS Collaboration, Search for ttH production using the Matrix Element Method, Tech. 

Rep. CMS-PAS-HIG-14-010, CERN, Geneva, 2014. 

[75] CMS Collaboration, V. Khachatryan et. al., Observation of the diphoton decay of the Higgs 
boson and measurement of its properties, Eur.Phys.J. C74 (2014), no. 10 3076, 

[arXiv: 1407.0558]. 

[76] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymptotic formulae for likelihood-based 
tests of new physics, Eur.Phys.J. C71 (2011) 1554, arXiv: 1007.1727 1 . 

[77] ATLAS Collaboration, Projections for measurements of Higgs boson cross sections, 
branching ratios and coupling parameters with the ATLAS detector at a HL-LHC, Tech. Rep. 
ATL-PHYS-PUB-2013-014, CERN, Geneva, Oct, 2013. 

[78] UKQCD Collaboration, G. Bali et. al., A Comprehensive lattice study of SU(3) glueballs, 
Phys.Lett. B309 (1993) 378-384, [hep-lat/9304012 . 

[79] W. Ochs, The Status of Glueballs, J.Phys. G40 (2013) 043001, (arXiv: 1301.5183]. 

[80] Particle Data Group Collaboration, J. Beringer et. al., Review of Particle Physics (RPP), 
Phys.Rev. D86 (2012) 010001. 

[81] ATLAS Collaboration, G. Aad et. al., Triggering on Long-Lived Neutral Particles in the 


62 


ATLAS Detector, Tech. Rep. ATL-PHYS-PUB-2009-082, ATL-COM-PHYS-2009-233, 2009. 

[82] CMS Collaboration, Search for Heavy Resonances Decaying to Long-Lived Neutral Particles 
in the Displaced Lepton Channel, Tech. Rep. CMS-PAS-EXO-11-004, CERN, Geneva, 2011. 

[83] CMS Collaboration, S. Chatrchyan et. al, Search for invisible decays of Higgs bosons in the 
vector boson fusion and associated ZH production modes, Eur.Phys.J. C74 (2014), no. 8 
2980, jarXiv: 1404.1344], 

[84] N. Arkani-Hamed and N. Weiner, LHC Signals for a SuperUnified Theory of Dark Matter, 
JHEP 0812 (2008) 104, jarXiv:0810.0714], 

[85] Gfitter Collaboration, M. Baak et. al., The global electroweak fit at NNLO and prospects for 
the LHC and ILC, Eur.Phys.J. C74 (2014), no. 9 3046, arXiv : 1407.3792 1 . 

[86] B. Holdorn, Two U(l)’s and Epsilon Charge Shifts, Phys.Lett. B166 (1986) 196. 

[87] S. Davidson, S. Hannestad, and G. Raffelt, Updated bounds on millicharged particles, JHEP 
0005 (2000) 003, [hep-ph/0001179 . 

[88] R. Essig, J. A. Jaros, W. Wester, P. H. Adrian, S. Andreas, et. al., Working Group Report: 
New Light Weakly Coupled Particles, arXiv: 1311.0029, 

[89] CMS Collaboration, S. Chatrchyan et. al., Search for Physics Beyond the Standard Model 
Using Multilepton Signatures in pp Collisions at y/s = 7 TeV, Phys.Lett. B704 (2011) 
411-433, [arXiv: 1106.0933]. 

[90] CMS Collaboration, S. Chatrchyan et. al., Search for electroweak production of charginos 
and neutralinos using leptonic final states in pp collisions at y/s = 7 TeV, JHEP 1211 
(2012) 147, arXiv: 1209.6620], 

[91] ATLAS Collaboration, G. Aad et. al, Search for New Phenomena in Events with Three 
Charged Leptons at y/s = 7 TeV with the ATLAS Detector, Phys.Rev. D87 (2013) 052002, 

[arXiv: 1211.6312], 

[92] ATLAS Collaboration, G. Aad et. al., Search for Anomalous Production of Prompt 
Like-Sign Lepton Pairs at y/s = 7 TeV with the ATLAS Detector, JHEP 1212 (2012) 007, 

[arXiv: 1210.4538]. 

[93] DO Collaboration, V. Abazov et. al, Search for dark photons from supersymmetric hidden 
valleys, Phys.Rev.Lett. 103 (2009) 081802, jarXiv:0905.1478 1 . 

[94] DO Collaboration, V. M. Abazov et. al, Search for events with leptonic jets and missing 
transverse energy in pp collisions at y/s = 1.96 TeV, Phys.Rev.Lett. 105 (2010) 211802, 


63 



arXiv:1008.3356. 

[95] CMS Collaboration, S. Chatrchyan et. ai, Search for Light Resonances Decaying into Pairs 
of Muons as a Signal of New Physics, JHEP 1107 (2011) 098, arXiv: 1106.2375j. 

[96] CMS Collaboration, S. Chatrchyan et. al, Search for a non-standard-model Higgs boson 
decaying to a pair of new light bosons in four-muon final states, Phys.Lett. B726 (2013) 
564-586, [arXiv: 1210.7619]. 

[97] ATLAS Collaboration, A Search for Lepton-Jets with Muons at ATLAS, Tech. Rep. 
ATLAS-CONF-2011-076, CERN, Geneva, May, 2011. 

[98] ATLAS Collaboration, G. Aad et. ai, A search for prompt lepton-jets in pp collisions at 
y/s = 7 TeV with the ATLAS detector, Phys.Lett. B719 (2013) 299-317, arXiv: 1212.5409 . 

[99] CMS Collaboration, V. Khachatryan et. al., Search for long-lived particles that decay into 
final states containing two electrons or two muons in proton-proton collisions at yfs = 8 
TeV, arXiv: 1411.6977. 


64 


